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Abstract. 

We present the results of our detailed pscudospectral direct numerical 
simulation (DNS) studies, with up to 1024'^ collocation points, of incompressible, 
magnetohydrodynamic (MHD) turbulence in three dimensions, without a mean 
magnetic field. Our study concentrates on the dependence of various statistical 
properties of both decaying and statistically steady MHD turbulence on the magnetic 
Prandtl number PrM over a large range, namely, 0.01 < PrM < 10. We obtain data for 
a wide variety of statistical measures such as probability distribution functions (PDFs) 
of moduli of the vorticity and current density, the energy dissipation rates, and velocity 
and magnetic-field increments, energy and other spectra, velocity and magnetic-field 
structure functions, which we use to characterise intermittency, isosurfaces of quantities 
such as the moduli of the vorticity and current, and joint PDFs such as those of fluid 
and magnetic dissipation rates. Our systematic study uncovers interesting results 
that have not been noted hitherto. In particular, we find a crossover from larger 
intermittency in the magnetic field than in the velocity fleld, at large PrM, to smaller 
intermittency in the magnetic field than in the velocity field, at low PrM. Furthermore, 
a comparison of our results for decaying MHD turbulence and its forced, statistically 
steady analogue suggests that we have strong universality in the sense that, for a fixed 
value of PrM, multiscaling exponent ratios agree, at least within our errorbars, for both 
decaying and statistically steady homogeneous, isotropic MHD turbulence. 



PACS numbers: 47.27.Gs,47.65.-Fa,05.45.-a 
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1. Introduction 

The hydrodynamics of conducting fluids is of great importance in many terrestrial 
and astrophysical phenomena. Examples include the generation of magnetic fields via 
dynamo action in the interstellar medium, stars, and planets PQlSlElSllSliniElEl 
El [ini E], and liquid-metal systems [El 131 [H US HEl [HI IE] that are studied in 
laboratories. The flows in such settings, which can be described at the simplest level by 
the equations of magnetohydrodynamics (MHD), are often turbulent [5]. The larger the 
kinetic and magnetic Reynolds numbers. Re = UL/u and Rcm = UL/rj, respectively, 
the more turbulent is the motion of the conducting fluid; here L and U are typical 
length and velocity scales in the flow, u is the kinematic viscosity, and r] is the magnetic 
diffusivity. The statistical characterization of turbulent MHD flows, which continues to 
pose challenges for experiments [I9], direct numerical simulations ^Q\, and theory [2T] . 
is even harder than its analogue in fluid turbulence because (a) we must control both Re 
and Rcm and (b) we must obtain the statistical properties of both velocity and magnetic 
fields. 

The kinematic viscosity u and the magnetic diffusivity t] can differ by several 
orders of magnitude, so the magnetic Prandtl number PrM = ReM/Re = I'/rj can 
vary over a large range. For example, PrM — 10~^ in the liquid-sodium system |15l 116] . 
PrM — 10~^ at the base of the Sun's convection zone [22] . and PrM ^ lO^'' in the 
interstellar medium [HI [20]. Furthermore, two dissipative scales play an important role 
in MHD; they are the Kolmogorov scale id (~ u^^^ at the level of Kolmogorov 1941 
(K41) phenomenology [231 Elj) and the magnetic-resistive scale (~ 77^/'^ in K41). A 
thorough study of the statistical properties of MHD turbulence must resolve both these 
dissipative scales. Given current computational resources, this is a daunting task at large 
Re especially if PrM is significantly different from unity. Thus, most direct numerical 
simulations (DNS) of MHD turbulence [251 ESI EH EHl EH ED] have been restricted to 
PrM — 1. Some DNS studies have now started moving away from the PrM — 1 regime 
especially in the context of the dynamo problem [STl [52] . 

Here we initiate a detailed DNS study of the statistical properties of incompressible, 
homogeneous, and isotropic MHD turbulence for a large range of the magnetic Prandtl 
number, namely, 0.01 < PrM < 10. There is no mean magnetic field in our DNS [33] ; 
and we restrict ourselves to Eulerian measurements [34j. Before we give the details 
of our DNS study, we highlight a few of our qualitative, principal results. Elements 
of some of our results, for the case PrM = 1 and for quantities such as energy 
spectra, exist in the MHD-turbulence literature as can be seen from the representative 
Refs. [51 [6l [251 [SHI EHl [371 I3H]- However, to the best of our knowledge, no study 
has attempted as detailed and systematic an investigation of the statistical properties 
of MHD turbulence as we present here, especially with a view to elucidating their 
dependence on PrM. Our study uncovers interesting trends that have not been noted 
hitherto. These emerge from our detailed characterisation of intermittency, via a 
variety of measures which include various probability distribution functions (PDFs) 
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such as those of the modulus of the vorticity and the energy dissipation rates, velocity 
and magnetic-field structure functions that can be used to characterise intermittency, 
isosurfaces of quantities such as the moduli of the vorticity and current, and joint PDFs 
such as those of fluid and magnetic dissipation rates. Earlier DNS studies [30| have 
suggested that intermittency, as characterised say by the multiscaling exponents for 
velocity- and magentic-field structure functions, is more intense for the magnetic field 
than for the velocity field when PrM = 1. Our study confirms this and suggests, 
in addition, that this result is reversed as we lower PrM- This crossover from larger 
intermittency in the magnetic field than in the velocity field, at large PrM, to smaller 
intermittency in the magnetic field than in the velocity field, at low Pr^i, shows up not 
only in the values of multiscaling exponent ratios, which we obtain from a detailed 
local-slope analysis of extended-self-similarity (ESS) plots [SHI HO] of one structure 
function against another, but also in the behaviours of tails of PDFs of dissipation 
rates, the moduli of vorticity and current density, and velocity and magnetic-field 
increments. Furthermore, a comparison of our results for decaying MHD turbulence 
and its forced, statistically steady analogue suggest that, at least given our conservative 
errors, the homogeneous, isotropic MHD turbulence that we study here displays strong 
universality [ITl 112] in the sense that multiscaling exponent ratios agree for both 
decaying and statistically steady cases. 

The remaining part of this paper is organised as follows. In Sec. [2] we describe the 



MHD equations, the details of the numerical schemes we use (Subsection 2.1), and the 



statistical measures we use to characterise MHD turbulence (Subsection 2.2). In Sec. |3] 



we present our results; these are described in the seven Subsections |3.H [3l2| |3.3[ |3^ |3.5[ 



3.6, and 3.7 that are devoted, respectively, to (a) a summary of well-known results from 
fluid turbulence that are relevant to our study, (b) the temporal evolution of quantities 
such as the energy and energy-dissipation rates, (c) energy, dissipation-rate, Elsasser- 
variable, and effective-pressure spectra, (d) various probability distribution functions 
(PDFs) that can be used, inter alia, to characterise the alignments of vectors such as 
the vorticity with, say, the eigenvectors of the rate-of-strain tensor, (e) velocity and 
magnetic-field structure functions that can be used to characterise intermittency, (f) 
isosurfaces of quantities such as the moduli of the vorticity and current, and (g) and 
joint PDFs such as those of fluid and magnetic dissipation rates. Section |4] contains a 
discussion of our results. 

2. MHD Equations 

The hydrodynamics of a conducting fluid is governed by the MHD equations [H 121 El IH 
ISl [7] , in which the Navier-Stokes equation for a fluid is coupled to the induction equation 
for the magnetic field: 

^ + (u ■ V)u = z/V^u - Vp + (b ■ V)b + A, (1) 
(9b 

— + (u ■ V)b = (b ■ V)u + r^V^b + (2) 
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here u, b, a; = V x u and j = V x b are, respectively, the velocity field, the magnetic 
field, the vorticity, and the current density at the point x and time t; v and r/ are the 
kinematic viscosity and the magnetic diffusivity, respectively, and the effective pressure 
is p = p + (fo^/Svr), where p is the pressure; /„ and are the external forces; while 
studying decaying MHD turbulence we set = = 0. The MHD equations can also 
be written in terms of the Elsasser variables = u±b [ZIES]- We restrict ourselves to 
low-Mach-number fiows so we use the incompressibility condition V ■ u(x, t) = 0; and 
we must, of course, impose V ■ b(x, t) = 0. By using the incompressibility condition, 
we can eliminate the effective pressure and obtain the velocity and magnetic fields via a 
pseudospectral method that we describe in the next Subsection. The effective pressure 
then follows from the solution of the Poisson equation 



2.1. Direct Numerical Simulation 

Our goal is to study the statistical properties of homogeneous and isotropic MHD 
turbulence so we use periodic boundary conditions and a standard pseudospectral 
method |13] with collocation points in a cubical simulation domain with sides of 
length L = 27r; thus, we evaluate spatial derivatives in Fourier space and local products 
of fields in real space. We use the 2/3 dealiasing method ^Sj to remove aliasing errors; 
after this dealiasing /cmax is the magnitude of the largest-magnitude wave vectors resolved 
in our DNS studies. We have carried out extensive simulations with = 512 and 
N = 1024; the parameters that we use for different runs are given in Table [T] for 
decaying and statistically steady turbulence. 

We use a second-order, slaved, Adams-Bashforth scheme, with a time step 6t, for 
the time evolution of the velocity and magnetic fields; this time step is chosen such that 
the Courant-Friedrichs-Lewy (CFL) condition is satisfied 

In our decaying-MHD-turbulence studies we have taken the initial (superscript 0) 
energy spectra E'^{k) and E^{k), for velocity and magnetic fields, respectively, to be the 
same; specifically, we have chosen 



where i?^, the initial amplitude, is chosen in such a way that we resolve both fiuid 
and magnetic dissipation scales and rj^, respectively: in all, except a few, of our 
runs k^ajjf^ > 1 and k^a.^rf^ > 1. The initial phases of the Fourier components 
of the velocity and magnetic fields are taken to be different and chosen such that 
they are distributed randomly and uniformly between and 2tt. In such studies, it 
is convenient to pick a reference time at which various statistical properties can be 
compared. One such reference time is the peak that occurs in a plot of the energy 
dissipation versus time; this reference time has been used in studies of decaying fluid 
turbulence |l5l HB], decaying fluid turbulence with polymer additives [HI SH], and 
decaying MHD turbulence ESI SH] • Such peaks are associated with the completion of 



V2p = V-[(b-V)b-(u-V)u]. 



(3) 



El{k) = El{k) = E%''exp{-2e), 



(4) 



Systematics of the Ftm dependence of homogeneous, isotropic MHD turbulence 5 

Table 1. List of parameters for our 16 DNS runs R1-R5, R3B-R5B, R1C-R4C, and 
R1D-R4D: is the number of collocation points in our simulation, v is the kinematic 
viscosity, PrM is magnetic Prandtl number, St is time step; Wrms, A, and Rca are the 
root-mean-sqare velocity, the Taylor microscale, and the Taylor-microscale Reynolds 
number, respectively. These are obtained at tc for our decaying-MHD-turbulence runs 
R1-R5, R3B-R5B, and R1C-R4C; and for statistically steady MHD turbulence (runs 
R1D-R4D) these are averaged over the statistically steady state; here tc (iteration steps 
multiplied by 5t) is the time at which the cascades for both the fluid and magnetic 
fields are completed (see text); rj^ and 77^ are, respectively, the Kolmogorov dissipation 
length scales for the fluid and magnetic fields, fcmax is the magnitude of the largest- 
magnitude wave vectors resolved in our DNS studies which use the 2/3 dealiasing rule; 
fcmax ^ 170.67 and 341.33 for N = 512 and 1024, respectively. 



Runs 




V 


PrM 


bt 


^rms 


A 


ReA 








Rl 


512 


10-4 


0.1 


10-3 


0.34 


0.18 


610 


9.7 


0.629 


2.280 


R2 


512 


5 X 10-4 


0.5 


10-3 


0.34 


0.27 


187 


9.1 


1.752 


2.389 


R3 


512 


10-3 


1 


10-3 


0.34 


0.35 


121 


8.1 


2.772 


2.444 


R4 


512 


5 X 10-3 


5 


10-3 


0.33 


0.60 


39 


7.0 


8.224 


2.692 


R5 


512 


10-2 


10 


10-3 


0.31 


0.73 


23 


6.5 


13.267 


2.836 


R3B 


512 


10-3 


1 


10-4 


1.07 


0.20 


210 


3.1 


1.175 


1.052 


R4B 


512 


5 X 10-3 


5 


10-4 


2.32 


0.24 


110 


1.4 


1.961 


0.644 


R5B 


512 


10-2 


10 


10-4 


3.21 


0.26 


85 


1.0 


2.490 


0.520 


RIC 


1024 


10-4 


0.01 


10-4 


0.35 


0.23 


810 


8.0 


1.431 


22.12 


R2C 


1024 


10-4 


0.1 


10-4 


1.11 


0.08 


890 


2.9 


0.472 


1.690 


R3C 


1024 


10-3 


1 


10-4 


1.14 


0.15 


172 


2.5 


1.996 


1.779 


R4C 


1024 


10-2 


10 


10-4 


2.37 


0.24 


57 


1.1 


5.550 


1.164 


RID 


512 


10-4 


0.01 


10-4 


1.31 


0.18 


2367 




0.320 


5.364 


R2D 


512 


10-4 


0.1 


10-4 


0.99 


0.14 


1457 




0.334 


1.145 


R3D 


512 


10-3 


1 


10-4 


1.06 


0.17 


239 




1.264 


1.033 


R4D 


512 


10-2 


10 


10-4 


1.04 


0.23 


61 




6.505 


1.129 



the energy cascade from large length scales, at which energy is injected into the system, 
to small length scales at which viscous losses are significant. In the MHD case, these 
peaks occur at slightly different times, t„ and tf,, respectively, in plots of the kinetic (e„) 
and magnetic (efe) energy-dissipation rates. In our decaying-MHD-turbulence studies 
we store velocity and magnetic fields at the time tc] if tu > tb, tc = tu', and tc = tf, 
otherwise; from these fields we calculate the statistical properties that we present in the 
next Section. 

In the simulations in which we force the MHD equations to obtain a nonequilibrium 
statistically steady state (NESS), we use a generalization of the constant-energy- 
injection method described in Ref. [50]. We do not force the magnetic field directly 
so we choose = 0. The force /„(x, t) is specified most simply in terms of /„(k, t), its 
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spatial Fourier components, as follows: 

where Q{kf — k) is 1 if k < kj and otherwise, V is the power input, and Eu{kf,t) = 
'^k<kf i); in our DNS we use kf = 2. This yields a statistically steady state in 
which the mean value of the total energy dissipation rate per unit volume balances the 
power input, i.e., 

(6) = V; (6) 

once this state has been established, we save 50 representative velocity- and magnetic- 
field configurations over ~ 36.08^/, 29.29t/, 32.61t/, and 30.95t/, for RID, R2D, R3D, 
and R4D, respectively, where tj = ii/urms is the integral-scale eddy-turnover time. We 
use these configurations to obtain the statistical properties that we describe below. 

For decaying MHD turbulence we have carried out eight simulations with 512^ 
collocation points and four simulations with 1024'^ collocation points. The parameters 
used in these simulations, which we have organised into three sets, are given in Table [T] 

In the first set of runs, R1-R5, we set the magnetic diffusivity rj = 10~^ and 
use five values of z/, namely, 10~^, 5.0 x 10~^, 10~^, 5.0 x 10"'^, and 10^^, which yield 
PrM = 0.1, 0.5, 1, 5, and 10. These runs have been designed to study the effects, on 
decaying MHD turbulence, of an increase in PrM, with the initial energy held fixed: in 
particular, we use = ~ 0.32 in Eq. |4] for runs R1-R5. Given that this initial 
energy and rj are both fixed, an increase in PrM leads to a decrease in Re, and thus an 
increase in k^.^jrf^ and k^^^^rf^ as we discuss in detail later. 

In our second set of decaying-MHD-turbulence runs, R3B, R4B and R5B, we 
increase E'^ in Eq. |4] as we increase z/, and thereby PrM, so that kj^g^^Vd — ^ 
kmaxVd — ^- Thus, in these runs, the inertial ranges in energy spectra extend over 
comparable ranges of the wave- vector magnitude k. 

Our third set of decaying-MHD-turbulence runs, RIC, R2C, R3C, and R4C, use 
1024^ collocation points and PrM = 0.01, 0.1, 1, and 10, respectively. By comparing the 
results of these runs with those of R1-R5, R3B, R4B, and R5B, we can check whether 
our qualitative results depend significantly on the number of collocation points that we 
use. 

We have carried out another set of four runs, RID, R2D, R3D, and R4D, in which 
we force the MHD equations, as described above, until we obtain a nonequilibrium 
statistically steady state. These runs help us to compare the statistical properties of 
decaying and statistically steady turbulence. In these runs we use 512^ collocation 
points, and u and rj such that PrM = 0.01, 0.1, 1, and 10, respectively. 



2.2. Statistical measures 

We use several statistical measures to characterise homogeneous, isotropic MHD 
turbulence. Some, but not all, of these have been used in earlier DNS studies [251 
ESI ESI EH EB ESI IS3] and solar- wind turbulence [Ml E51 ES] • 
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We calculate the kinetic, magnetic, and total energy spectra Eu{k) = 
Ek9|k|=fc|u(k)P, Ebik) = Ek9|k|=fc|b(k)|2, and E{k) = E^ik) + Eb{k), respectively, 
the kinetic, magnetic, and total energies Eu = J2k Eu{k)/2, E^ = J2k-Eb{k)/2, and 
E = Eu + Eb, and the ratio Eb/Eu. Spectra for the Elsasser variables, energy dissipation 
rates, and the effective pressure are, respectively, Ez±{k) = X]k9|k|=fc I^^C^)!^; ^u{k) = 
uk'E^ik), eb{k) = uk'Ebik), and P{k) = Ek9|k|=fc 1^1'- 

Our MHD simulations are characterised by the Taylor-microscale Reynolds number 
Rca = MrmsA/i^, the magnetic Taylor-microscale Reynolds number RmA = MrmsA/77, 
and the magnetic Prandtl number PrM = Rm^/Re^ = i^/tj, where the root-mean- 
square velocity Uj-ms = ■\/2i?„/3 and the Taylor microscale A = [J2k k'^E{k)/E] ^^"^ . We 
also calculate the integral length scale ii = Efc-^(^)/^] /-^' mean kinetic energy 
dissipation rate per unit mass, = J^'^ijidiUj + djUiY = uj^k^'^^uik), the mean 
magnetic energy dissipation rate per unit mass tb = V J2i ji^i^j = V J2k k'^^bik), 

the mean energy dissipation rate per unit mass e = + tb, and the dissipation length 
scales for velocity and magnetic fields rj^ = (i^^/e„)^''^ and r/^ = {rj^/ebY^'^, respectively. 

We calculate the eigenvalues A'^ and the associated eigenvectors e"^, with n = 1,2, 
or 3, of the rate-of-strain tensor § whose components are Sij = diUj + djUi. Similarly 
A\, A2 and A| denote the eigenvalues of the tensile magnetic stress tensor T, which has 
components Tij = —bibj; the corresponding eigenvectors are, respectively, e\, 63, and Cg. 

For incompressible flows En^n = O5 so at least one of the eigenvalues A"^ must 
be positive and another negative; we label them in such a way that A3 is positive. A" 
negative, and Ag lies in between them; note that A2 can be positive or negative. We 
obtain probability distribution functions (PDFs) of these eigenvalues; furthermore, we 
obtain PDFs of the cosines of the angles that the associated eigenvectors make with 
the vectors such as u, u:, etc. These PDFs and those of quantities such as the local 
cross helicity He = u ■ b help us to quantify the degree of alignment of pairs of vectors 
such as u and b [52]. We also compare PDFs of magnitudes of local vorticity u, the 
current density j, and local energy dissipation rates and to get information about 
intermittency in velocity and magnetic fields. 

We also obtain several interesting joint PDFs; to the best of our knowledge, these 
have not been obtained earlier for MHD turbulence. We first obtain the velocity- 
derivative tensor A, also known as the rate-of-deformation tensor, with components 
^ij = diUj, and thence the invariants Q = — ^tr(A^) and R = — |tr(A^), which have 
been used frequently to characterise fluid turbulence [571 EHl EH] . The zero-discriminant 
line D = f i?^ + q3 ^ q ^^le Q and R axes divide the QR plane into qualitatively 
different regimes. In particular, regions in a turbulent flow can be classified as follows: 
when Q is large and negative, local strains are high and vortex formation is not favoured; 
furthermore, if i? > 0, fluid elements experience axial strain, whereas, if i? < 0, they feel 
biaxial strain. In contrast, when Q is large and positive, vorticity dominates the flow; 
if, in addition, R < 0, vortices are compressed, whereas, if R > 0, they are stretched. 
Thus, some properties of a turbulent flow can be highlighted by making contour plots 
of the joint PDF of Q and R; these QR plots show a characteristic, tear-drop shape. 
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We explore the forms of these and other joint PDFs, such as joint PDFs of eu and e^, 
in MHD turbulence. 

To characterise intermittency in MHD turbulence we calculate the order-p, equal- 
time, longitudinal structure functions Sp{l) = {\6a\\{'x, , where the longitudinal 
component of the field a is given by 6a^\ (x, /) = [a(x + 1, t) — a(x, t)] ■ |, where a can be 
u, b, or one of the Elsasser variables. From these structure functions we also obtain the 
hyperflatnessFg"(/) = S^{l)/[S^{l)f. For separat ions I in the inertial range, i.e., rj^, rj^ <^ 
/ ^ L, we expect Sp{l) ~ I'^p , where (p are the inertial-range multiscaling exponents 
for the field a; the Kolmogorov phenomenology of 1941 [23l [2^ [25] , henceforth referred 
to as K41, yields the simple scaling result Cp^"^^ = but multiscaling corrections are 
significant with 7^ Cp^^^ [Sec.jsj. From the increments ^ay (x, /) = [a(x+l, t)— a(x, t)]- 1 
we also obtain the dependence of PDFs of 6a\\ on the scale I. 

3. Results 

To set the stage for the types of studies we carry out for MHD turbulence, we 
begin with a very brief summary of similar and well-known results from studies 
of homogeneous, isotropic Navier-Stokes turbulence, which can be found, e.g., in 

Refs. [SassiliniiolEIlEaEsllSlESlE^ 

3.1. Overview of fluid turbulence 

For ready reference we show here illustrative plots from a DNS study that we have 
carried out for the three-dimensional Navier-Stokes equation by using a pseudospectral 
method, with 512^ collocation points and the 2/3 rule for removing aliasing errors; here 
1/ = 0.001, Rca ~ 340 and k^^^ri"^ ~ 0.3. 




Figure 1. Plots from our DNS of decaying fluid turbulence in the Navier-Stokes 
equation with 512'^ collocation points: (a) plots of the energy E (red full line) and 
mean energy dissipation rate e (blue dotted line) versus time t (given as a product 
of the number of iterations and the time step 5t); (b) log-log (base 10) plots of the 
energy spectrum E{k) (red dashed line) and the corresponding compensated spectrum 
k^/^Eik) (blue dotted line) versus k. The black solid line shows the K41 result k ^'^ 
for comparison; (c) log-log (base 10) plot of spectrum of energy-dissipation or enstrophy 
spectrum e(fc); (d) log-log (base 10) plots of the pressure spectrum P{k) (red dashed 
line) and the compensated pressure spectrum k'^/^P{k) (blue dotted line). The black 
solid line shows the K41 result k^'^/'^ for comparison. 
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Figure 2. PDFs from our DNS of decaying fluid turbulence in Navier-Stokes equation 
with 512'^ collocation points: (a) semilog (base 10) plots of PDFs of eigenvalues of 
rate-of-strain tensor §, namely, A" (red full line), (green dashed line) and A3 (blue 
dotted line); (b) semilog (base 10) plots of PDFs of cosines of angles between the 
vorticity a; and eigenvectors of S, namely, e" (red full line), 62 (green dashed line), and 
63 (blue dotted Hue); (c) semilog (base 10) plots of PDFs of cosines of angles between 
the velocity u and eigenvectors of S, namely, e" (red full line), 62 (green dashed line), 
and 63 (blue dotted line); (d) semilog (base 10) plots of PDFs of cosines of angles 
between the velocity u and vorticity u). 




Figure 3. PDFs from our DNS of decaying fluid turbulence in Navier-Stokes equation 
with 512"^ collocation points: semilog (base 10) plots of the PDFs of (a) the pressure 
p, (b) modulus of the vorticity uj and (c) the local energy-dissipation rate e. 




Figure 4. Isosurfaces of (a) the modulus of the vorticity w, (b) the local energy- 
dissipation rate e, and (c) the local pressure p, from our DNS of decaying fluid 
turbulence in Navier-Stokes equation with 512'^ collocation points. The isovalues used 
in these plots are two standard deviations more than the mean values of the quantities. 



In decaying fluid turbulence, energy is injected at large spatial scales as described 
in the previous Section for the MHD case. This energy cascades down till it reaches 
the dissipative scale at which viscous losses are significant. We study various statistical 
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Figure 5. QR plot, i.e., the joint PDF of Q and R (see text) shown as filled contour 
plot in our log-log (base 10) scale, obtained from a DNS of decaying fluid turbulence 
in Navier-Stokes equation with 512'^ collocation points. 



properties; these are given in points (i)-(vi) below: 

(i) Plots of the energy E and the mean energy dissipation rate e versus time show, 
respectively, a gentle decay and a maximum as shown, e.g., by the full red and dotted 
blue curves in Fig. |3l|(a). This maximum in e is associated with the completion of 
the energy cascade at a time tc, the remaining properties (ii)-(vi) are obtained at t^. 
(ii) If Rca is sufficiently large and we have a well- resolved DNS (i.e., /cmax^^ > 1); 
then, at tc, the spectrum E{k) shows a well-developed inertial range, where at the 
K41 level E{k) ~ k~^^^, and a dissipation range, in which the behaviour of the energy 
spectrum is consistent with E{k) ~ fc"exp(— /3fc), where a and P are non-universal, 
positive constants [EH El] and 5kd < k < lOkd, with ka = An illustrative energy 

spectrum is shown by the dashed red line in Fig. 3.1[ b); the blue dotted curve shows 
the compensated spectrum k^^^E{k); the associated dissipation or enstrophy spectrum 
e{k) is shown in Fig. 3.1 c) and the inertial-range pressure spectrum [67J, P{k) ~ /c"^/^ 



at the K41 level is shown in Fig. 3.1[ d). [Note that our DNS for the Navier-Stokes 
equation, which suffices for our purposes of illustration, does not have a well-resolved 
dissipation range because kma^Vd — 0-3 < 1; this is also reflected in the lack of a 
well-developed maximum in the enstrophy spectrum of Fig. 3.1 c).] (iii) Illustrative 



PDFs of the eigenvalues A" of the rate-of-strain tensor § are given for n = 1, 2, and 
3, respectively, by the full red, dashed green, and dotted blue curves in Fig. |3.lK a); 
PDFs of the cosines of the angles that the vorticity and the velocity u make with 
the associated eigenvectors ej^ are given, respectively, in Figs. 3.1[ b) and 3.1[ c) via full 
red (ra = 1), dashed green {n = 2), and dotted blue {n = 3) curves; these show that 
both u} and u have a tendency to be preferentially aligned parallel or antiparallel to 
^2 [Sn]; the PDF of the cosine of the angle between u and also indicates preferential 
alignment or antialignment of these two vectors, but with a greater tendency towards 
alignment as found in experiments with a small amount of helicity |68j and as illustrated 
m our Fig. |3ljd). Finally, we give representative PDFs of the pressure p, modulus of 
vorticity cu = and the local energy dissipation e in Figs. 3.1[a), 3.1[ b), and 3.1[ c), 
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respectively; note that the PDF of the pressure is negatively skewed, (iv) Inertial- 
range structure functions Sp{l) ~ l^p show significant deviations [21] from the K41 
result Cp^'^^ = p/3 especially for p > 3. From these structure functions we can obtain 
the hyperflatness Fq{1); this increases as the length scale / decreases, a clear signature 
of intermittency, as shown, e.g., in Refs. |1HI E5]- This intermittency also leads to 
non-Gaussian tails, especially for small /, in PDFs of velocity increments (see, e.g., 
Refs. [HSl ESI EO]) such as Su\\{l). (v) Small-scale structures in turbulent flows can be 
visualised via isosurfaces [71] of, say, co, e, and p, illustrative plots of which are given 
in Figs. 3.1[ a)- |3A| ^c); these show that regions of large u are organised into slender 
tubes whereas isosurfaces of e look like shredded sheets; pressure isosurfaces also show 
tubes [361 SH] but some studies have suggested the term cloud-like for them [60]. (vi) 
Joint PDFs also provide useful information about turbulent flows; in particular, contour 
plots of the joint PDF of Q and R, as in the representative Fig. |5} show a characteristic 
tear-drop structure. 

The properties of statistically steady, homogeneous, isotropic fluid turbulence are 
similar to those described in points (ii)-(vi) in the preceding paragraph for the case 
of decaying fluid turbulence at cascade completion at t^. In particular, the strong- 
universality |41| hypothesis suggests that the multiscaling exponents have the same 
values in decaying and statistically steady turbulence. 

The remaining part of this Section is devoted to our detailed study of the MHD- 
turbulence analogues of the properties (i)-(vi) summarised above; these are discussed, 
respectively, in the six Subsections 3.2 3.7[ 



3.2. Temporal evolution 

We examine the time evolution of the energy, the energy-dissipation rates, and related 
quantities, first for decaying and then for statistically steady MHD turbulence. 

Figure [6] shows how the total energy E (red full line), the kinetic energy Ey_ (green 
dashed line), and the magnetic energy (blue dotted line) evolve with time t (given as 
a product of the number of iterations and the time step St) for runs R1-R5 [Figs. [6]^a.l)- 
|6]^e.l)] and runs R3B-R5B [Figs. [6]^f.l)-|6](h.l)] for decaying MHD turbulence. Figure |6] 
also shows similar plots for the mean energy dissipation rate e (red full line), the mean 
kinetic-energy dissipation rate (green dashed line), and the mean magnetic-energy 
dissipation rate eb (blue dotted line) versus time t for runs R1-R5 [Figs.|6|a.2)-[6|e.2)] and 
runs R3B-R5B [Figs. [6]^f.2)-|6](h.2)]. In addition Fig. |6] depicts the time-evolution of the 
ratio Eb/E^ for runs R1-R5 [Figs. ga.3)-ge.3)] and runs R3B-R5B [Figs. gf.3)@h.3)]. 
We see from these figures that, for all the values of PrM we have used, the energies E 
and Eu decay gently with t but Eb rises initially such that the ratio Eb/E^ rises, nearly 
monotonically, with t over the times we have considered; this is an intriguing trend that 
does not seem to have been noticed earlier. The times over which we have carried out 
our DNS are comparable to the cascade-completion time tc that can be obtained from 
the peaks in the plots of e, e„, and versus t [Figs. |6](a.2)-|6](h.2)]; by comparing these 



Systematics of the Piu dependence of homogeneous, isotropic MHD turbulence 



12 



(c.l)- 



0.7 






0.6 






0.5 






0.4 






0.3 






0.2 












' 0.0 1.5 

70 


3.0 4.5 6.0 7.5 9.0 


60 
50 




(h.l); 


40 
30 






20 
10 






n 








Figure 6. Plots versus time t (given as a product of tlic number of iterations 
and the time step St) of energies (a.l)-(h.l); total energy E (red full line), kinetic- 
energy Eu (green dashed line) , and magnetic-energy Ef, (blue dotted line) , of energy- 
dissipation rates (a.2)-(h.2); mean energy dissipation rate e (red, full line), kinetic- 
energy dissipation e„ (green dashed line), and magnetic-energy dissipation rate (blue 
dotted line), and of the ratio E^/Eu (a.3)-(h.3), generically, for decaying simulations 
(a) PrM = 0.1 (Rl), (b) PrM = 0.5 (R2), (c) Prw = 1.0 (R3), (d) Prw = 5.0 (R4), 
(e) PrM = 10.0 (R5), (f) Prw = 1.0 (R3B), (g) Prw = 5.0 (R4B), and (h) PrM = 10.0 
(R5B). 



plots we see that, as we move from Pim = 0.1 to Pim = 10, with fixed rj, we find that 
(e„ — ef,) and (t;, — t„) grow from negative values to positive ones because e„ increases 
with PrM, where th and tu are the positions of the cascade-completion maxima in eh 
and e^i respectively. We do not pursue the time evolution of our system well beyond t^ 
and tb because the integral scale begins to grow thereafter and, eventually, can become 
comparable to the linear size of the simulation domain |45j . 

Figures [7]^a.l)-[7||^d.l), show how the total energy E (red full line), the total kinetic 
energy (green dashed line), and the total magnetic energy E^ (blue dotted line) 
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Figure 7. Plots versus time t (given as a product of tlic number of iterations 
and the time step St) of energies (a.l)-(d.l); total energy E (red full line), kinetic- 
energy Eu (green dashed line), and magnetic-energy Eb (blue dotted line), of energy- 
dissipation rates (a.2)-(d.2); mean energy dissipation rate e (red, full line), kinetic- 
energy dissipation rate e„ (green dashed line), and magnetic-energy dissipation rate Cb 
(blue dotted line), and of the ratio Eb/E^ (a.3)-(d.3), generically, for forced simulations 
(a) PrM = 0.01 (RID), (b) PrM = 0.1 (R2D), (c) Ptm = 1-0 (R3D), and (d) PrM = 10 
(R4D). 



evolve with time t (given as a product of the number of iterations and the time step 
6t) for, respectively, runs R1D-R4D for forced and statistically steady MHD turbulence. 
Figures [7](a.2)-[7](d. 2), show similar plots for the mean energy dissipation rate e (red full 
line), the mean kinetic-energy dissipation rate eu (green dashed line), and the mean 
magnetic-energy dissipation rate eb (blue dotted line) versus time t for, respectively, 
runs R1D-R4D. And Figs. [7](a.3)-[7](d.3), depict the time-evolution of the ratio E^/Eu 
for these runs. We see from these figures that a statistically steady state is established 
in which the energies E, Eu, and £^6, the dissipation rates e, e„, and e?,, and the ratio 
Eh/Eu fluctuate about their mean values (after initial transients have died out). The 
mean value of Eb/Eu increases from about 0.2 — 0.3 to Eb/Eu ~ 1 as PrM increases from 
0.01 to 10. Furthermore, the mean values of the dissipation rates and e^, are such 
that (e„ — Cfc) grows from a negative value ^ —1 to a value close to zero as PrM increases 
from 0.01 to 10. 
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Figure 8. Log-log (base 10) plots of the compensated energy spectra e "^/^k^/^E^k) 
(red full lines), e^2/3k^^^K{k) ( green dashed lines), and e '^/^k^/^Ei,{k) (blue dotted 
lines); on the vertical axes these are denoted generically as E'^{k): (a.l) PrM = 0.1 
(Rl), (b.l) PrM = 0.5 (R2), (c.l) Prw = 1.0 (R3), (d.l) Prw = 5.0 (R4), (e.l) 
PrM = 10.0 (R5), (f.l) PrM = 1.0 (R3B), (g.l) PrM = 5.0 (R4B), (h.l) Prw = 10.0 
(R5B), (a.2) PrM = 0.01 (RIC), (b.2) PrM = 0.1 (R2C), (c.2) PrM = 1.0 (R3C), and 
(d.2) PrM = 10.0 (R4C) for decaying MHD turbulence; and for statistically steady 
MHD turbulence (a.3) PrM = 0.01 (RID), (b.3) PrM = 0.1 (R2D), (c.3) Prw = 1.0 
(R3D), and (d.3) PrM = 10.0 (R4D). 



3.3. Spectra 

We now discuss the behaviours of the energy, kinetic-energy, magnetic-energy, Elsasser 
variable, dissipation-rate, and effective-pressure spectra, first for decaying and then for 
statistically steady MHD turbulence. In the former case, spectra are obtained at the 
cascade-completion time tc] in the latter, they are averaged over the statistically steady 
state that we obtain. 

We present compensated spectra of the total energy Ec{k) = e~'^^^k^^^E(k) (red 
full line), the kinetic energy E^{k) = e~'^^^k^^^Eu{k) (green dashed line), and the 
total magnetic energy El{k) = e-^/^k^/^Ek{k) (blue dotted line) at t^ for runs R1-R5 
[Figs. ga.l)-ge.l)], R3B-R5B [Figs. |8|(f.l)-gh.l)], and R1C-R4C [Figs. ga.2)-gd.2)] 
for decaying MHD turbulence; and runs R1D-R4D [Figs. |8]^a.3)-|8](d.3)] show these for 
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Figure 9. Log-log (base 10) plots of compensated energy spectra, Ef (k) = k^^^E^{k), 
with k the magnitude of the wave vector, for the Elsasser variables fields z+ (red full 
line) and (blue dashed line): (a.l) Ptm = 0.1 (Rl), (b.l) Pym = 0.5 (R2), (c.l) 
PrM = 1.0 (R3), (d.l) PrM = 5.0 (R4), (c.l) Prw = 10.0 (R5), (f.l) Prw = 1.0 (R3B), 
(g.l) PrM = 5.0 (R4B), (h.l) Prw = 10.0 (R5B), (a.2) PrM = 0.01 (RIC), (b.2) 
PrM = 0.1 (R2C), (c.2) Ptm = 1.0 (R3C), and (d.2) Ptm = 10.0 (R4C) for decaying 
MHD turbulence; and for statistically steady MHD turbulence (a.3) PrM = 0.01 (RID), 
(b.3) PrM = 0.1 (R2D), (c.3) Ptm = 1.0 (R3D), and (d.3) PrM = 10.0 (R4D). 

statistically steady MHD turbulence. From Figs. [8]^a.l)-[8]^e.l) and Table [l] we see that 
increases as we increase u to increase PrM, because the initial energy is the same for 
runs R1-R5, so the dissipation tail in E^{k) is drawn in towards smaller and smaller 
values of k as we move from PrM = 0.1 to PrM = 10; between PrM = 0.5 and PrM = 1 the 
tails of E'^{k) and E^ik) and eventually E^{k) dominates and becomes indistinguishable 
from Ec{k) on the scales of Figs. |8](d.l) andjsj^e.l). A comparison of Figs. [8]^f.l)-|8](h.l) 
shows that, if we increase PrM from 1 to 10, we can keep both kmn^Vd ^max''?^ close 
to 1 so the dissipation ranges of these spectra span comparable ranges of k; however, 
as PrM increases, more and more of the energy is concentrated in the magnetic field. 
These trends are not affected (a) if we increase the number of collocation points, as can 
be seen from the compensated spectra in Figs. [8|^a.2)-[8]^d.2) for runs R1C-R4C, which 
use 1024^ collocation points, or (b) if we study energy spectra for statistically steady 



Systematics of the Ftm dependence of homogeneous, isotropic MHD turbulence 16 




Figure 10. Log-log (base 10) plots of energy-dissipation spectra for the fluid (red 
full lines) and magnetic (blue dashed lines) fields, with k the magnitude of the wave 
vector: (a.l) Ptm = 0.1 (Rl), (b.l) Ptm = 0.5 (R2), (c.l) Ptm = 1-0 (R3), (d.l) 
PrM = 5.0 (R4), (e.l) PrM = 10.0 (R5), (f.l) PrM = 1.0 (R3B), (g.l) Pim = 5.0 
(R4B), (h.l) PrM = 10.0 (R5B), (a.2) Pim = 0.01 (RIC), (b.2) Prw = 0.1 (R2C), (c.2) 
PrM = 1.0 (R3C), and (d.2) PrM = 10.0 (R4C) for decaying MHD turbulence; and for 
statistically steady MHD turbulence (a.3) PrM = 0.01 (RID), (b.3) PrM = 0.1 (R2D), 
(c.3) PrM = 1.0 (R3D), and (d.3) Ptm = 10.0 (R4D). 



MHD turbulence as can be seen from the compensated spectra in Figs. |8j(a.3)-|8j(d.3) 
for runs R1D-R4D. Figures gel), gg.l), |8|(c.2), andgc.3), for runs R3 (Rca = 121), 
R3B (ReA = 210), R3C (ReA = 172), and R3D (ReA = 239), respectively, all lie in one 
column and all have PrM = 1; so they provide a convenient way of comparing the ReA 
dependence of these spectra with a fixed value of Pym = 1- All the spectra in subfigures 
of Fig. [sjhave been compensated by the 5/3 power of k and, to the extent that they show 
small, flat parts, their inertial-range, energy-spectra scalings are consistent with k~^^^ 
behaviours; other powers, such as —3/2, can also give small, fiat parts in compensated 
spectra. A detailed error analysis is required to decide which power is most consistent 
with our data; we defer such an error analysis to the Subsection 3^ where we carry it 
out for structure functions. 

Compensated spectra of the Elsasser variables, namely, E^(k) = e^'^^^k^^^E~^(k) 
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Figure 11. (a) The kinetic energy spectrum Eu(k) (red asterisks) deep in the 
dissipation range for run R5; the black hne indicates the fit Eu{k) ~ k"^-^^ exp(— 0.235 k) 
for 5A;^ < k < lOfc^, where fc^ = l/??^; (b) the magnetic energy spectrum £'b(fc) (red 
asterisks) deep in the dissipation range for run RIC; the black line indicates the fit 
Eb{k) - fc-5-24exp(-0.014 k) for bk\ < k < lOfc^J, where fc^J = l/rf'^. 




Figure 12. Log-log (base 10) plots of effective pressure spectra P{k) (red full lines), 
with k the magnitude of the wave vector, and the corresponding compensated spectra 
P{k)ky^ (blue, dashed lines): (a.l) Ptm ^ 0.1 (Rl), (b.l) Ptm = 0.5 (R2), (c.l) 
PrM = 1.0 (R3), (d.l) PrM = 5.0 (R4), (e.l) Ptm = 10.0 (R5), (f.l) Ptm = 1-0 (R3B), 
(g.l) PrM = 5.0 (R4B), (h.l) Ptm = 10.0 (R5B), for decaying MHD turbulence; and 
for statistically steady MHD turbulence (a.2) PrM = 0.01 (RID), (b.2) PrM = 0.1 
(R2D), (c.2) PrM = 1.0 (R3D), and (d.2) PrM = 10.0 (R4D). 

(red full lines) and E~{k) = eu'^^^k^^^E~ (k) (blue dashed lines) are shown, at 
the cascade-completion time tc, for the decaying-MHD-turbulence runs R1-R5 in 
Figs.|9|(a.l)-|9|e.l), R3B-R5B in Figs. gf.l)-|9|h.l), and R1C-R4C in Figs. ga.2)-gd.2); 
and Figs. |8](a.3)-[9td.3) show these spectra for statistically steady MHD turbulence in 
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runs R1D-R4D, respectively. Note that the dissipation ranges of E^{k) and E~{k) 
overlap nearly on the scales of these figures. Differences between these are most 
pronounced at small /c, where, typically, E~{k) lies below E^{k); these differences 
decrease with increasing Prn, if we hold the initial energy fixed as in Figs. [9]^a.l)-[9]^e.l) 
for runs R1-R5. 

Next we come to the energy-dissipation (or enstrophy) spectra eu{k) = k^Eu{k) (red 
full line) and eb{k) = k'^E^^k) (blue dashed line) at tc- These are shown, at the cascade- 
completion time tc, for the decaying-MHD-turbulence runs R1-R5 in Figs. IToFa.l)- 



lOe.l), R3B-R5B in Figs. 10 ;f.l)-10;h.l), and R1C-R4C in Figs. Ma.2)-nord.2); and 



Figs. [To|^a.3)-|To|(d.3) depict these spectra for statistically steady MHD-turbulence runs 
R1D-R4D. To the extent that most of these spectra show maxima at values of k at 
the beginning of the dissipation range, most of our runs have well-resolved dissipation 
ranges; this also follows from the values of fcmax''?^ and kj^a^Vd Table [l} Runs RID and 
R2D have slightly under-resolved fluid-dissipation ranges with kma^Vd — 0-32 and 0.33, 
respectively; and, for the former, a barely discernible, dissipation-range maximum in 
e„(fc); however, as shown in our Navier-Stokes DNS in Subsection 3.1 , reasonable results 
can be obtained for various statistical properties with fcmax^d — 0.3. The elucidation 
of the behaviours of dissipation-range spectra of course require large values of fcmax^rf 
or k^gjjf^] in particular, runs R5 and RIC, with /cmax''?^ — 13.3 and k^gj^rf^ ~ 22.1, 
respectively, are well suited for uncovering the functional forms of Eu{k) and Eb{k) 
in their dissipation ranges. In Figs. |3.3[ a) and |3.3K b) we show, respectively, the 
kinetic- and magnetic-energy spectra Euik) and Ehik) deep in their dissipation ranges 
for runs R5 and Rl, respectively; our data for these spectra can be fit to the form 
~ A;" exp(— /3A;) for k deep in the dissipation range and a and /3 nonuniversal numbers 
that depend on the parameters of the simulation; similar results have been obtained for 
fluid turbulence [HIl El]- In particular, our data [Figs. |3.3K a) and 3.3[ b)] for runs R5 
and RIC are consistent with Eu{k) ~ A;^-''^ exp(-0.235 k), for ^k'^ < k < IQk^ with 
A;^ = l/ry^, and E^i^k) ~ A;-5-24 exp(-0.014 k), for 5k''^ < k < lOA;^ with A;^ = 1/r/^, 
respectively. 

We now turn to the spectra for the effective pressure P{k) (red full lines) and 
their compensated versions k'^/^P{k) (blue dashed lines) that are shown at t^ for runs 
R1-R5 [Figs. |l2ja.l)-|l2je.l)] and R3B-R5B [Figs. |l2jf.l)-|l2jh.l)] for decaying MHD 
turbulence; and for statistically steady MHD turbulence they are shown in Figs. [I2|^a.2)- 
[l2|d.2) for runs R1D-R4D. Pressure spectra have been studied for fluid turbulence as, 
e.g., in Refs. |16| [67]: to the best of our knowledge they have not been obtained for MHD 
turbulence. The compensated spectra here show that, for all our runs, the inertial- 
range behaviours of these effective-pressure spectra are consistent with the power law 
A;""^/^; this is consistent with the k~^/'^ behaviours of the energy spectra discussed above. 
Furthermore, as PrM increases from 0.1 to 10 in runs R1-R5, P{k) falls more and more 
rapidly as can be seen from the vertical scales in Figs. 12 a.l)-[l2|e.l). 
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3.4- Probability distribution functions 

We calculate several PDFs to characterise the statistical properties of decaying and 
statistically steady MHD turbulence. In the former case, PDFs are obtained at the 
cascade-completion time t^, in the latter, they are averaged over the statistically steady 
state that we obtain. The PDFs we consider are of two types: the first type are PDFs of 
the cosines of angles between various vectors, such as u and these help us to quantify 
the degrees of alignment between such vectors; the second type are PDFs of quantities 
such as eu, ^b, and the eigenvalues of the rate-of-strain tensor. 



In Fig. 13 we show plots of the PDFs of cosines of the angles between the vorticity 
iv and the eigenvectors of the fluid rate-of-strain tensor §, namely, (red full line), e„ 
(green dashed hues), and (blue dotted lines) for runs R1-R5 and R3B-R5B at the 
cascade-completion time tc for the case of decaying MHD turbulence. In Fig. [14] we 
show similar plots of the PDFs of cosines of the angles between the current density j 
and the eigenvectors of the fluid rate-of-strain tensor S. The most important features 
of these figures are sharp peaks in the green dashed lines; these show that there is a 
marked tendency for the alignment or antialignment of a; and e^, as in fluid turbulence, 
and of a similar tendency for the alignment or antialignment of j and e^; these features 
do not depend very sensitively on PrM. Furthermore, the PDFs of cosines of the angles 
between uj and (blue dotted lines) and uj and (red full lines) show peaks near zero 



in Fig. 13 in contrast, analogous PDFs for the cosines of the angles between j and 
(red full lines) and a; and (blue dotted lines) show nearly flat plateaux in the middle 



with very gentle maxima near —0.5 and 0.5 [Fig. 14 . Runs R1C-R4C and R1D-R4D 
yield similar PDFs, for the cosines of these angles, so we do not give them here. 

Plots of the PDFs of cosines of the angles between the velocity u and the 



eigenvectors of the fluid rate-of-strain tensor S are given Fig. 15; their analogues for 



b are given Fig. 16 Again, the most prominent features of these figures are sharp 
peaks in the green dashed lines; these show that there is a marked tendency for the 
alignment or antialignment of u and and of a similar tendency for the alignment or 
antialignment of j and e^; these features do not depend very sensitively on PrM. The 
PDFs of cosines of the angles between u and (red full line) and u and (blue dotted 
lines) show gentle, broad peaks that imply a weak preference for angles close to 45° 



or 135°; these peaks are suppressed as we increase Ftm [Figs. 15 a.l)-15 'e.l) for runs 



R1-R5] with fixed initial energy, but they reappear if we compensate for the increase of 



PrM by increasing the initial energy [Figs. 15 f.l)- 15 h.l)]. Similar, but sharper, peaks 



appear in the PDFs of cosines of the angles between u and (red full lines) and u 
and (blue dotted hues); these show a weak preference for angles close to 47° or 133° 
[Fig. [16]. Some simulations of compressible MHD turbulence have noted the presence 
of such peaks for PrM = 1. 

Only one of the eigenvalues of the tensile magnetic stress tensor T is non- 
zero; and the corresponding eigenvector el is identically aligned with b. Thus PDFs of 
cosines of angles between u, cj, j, and b and the eigenvectors of T are simpler than their 
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counterparts for § and are not presented here. 



Figure 17 shows plots of PDFs of cosines of angles, denoted generically by 9, between 

(a) u and b, (b) u and lo, (c) u and j, (d) uj and j, (e) b and u}, and (f) b and j for 
runs Rl (red lines), R2 (green lines), R3 (blue lines), R4 (black lines), and R5 (cyan 
lines). These figures show the following: (a) u and b are more aligned than antialigned 
[this is related to the small, positive, mean values of He (see below) in our runs R1-R5]; 

(b) u and u} and more antialigned than aligned, as noted for decaying fluid turbulence 
with slight helicity in Refs. |l6l EH]; (c) u and j show approximately equal tendencies 
for alignment and antialignment; (d) u: and j display a greater tendency for alignment 
than antialignment; (e) b and oj have approximately equal tendencies for alignment and 
antialignment; and (f) b and j are more antialigned than aligned. 
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Figure 13. Semilog (base 10) plots of the PDFs of cosines of the angles, denoted 
generically by 0, between the vorticity uj and the eigenvectors of the fluid rate-of-strain 
tensor S, namely, e„ (red full line),e^ (green dashed line), and (blue dotted line): 
(a.l) PrM = 0.1 (Rl), (b.l) Prw = 0.5 (R2), (c.l) Prw = 1.0 (R3), (d.l) Prw = 5.0 
(R4), (e.l) PrM = 10.0 (R5), (f.l) Prw = 1.0 (R3B), (g.l) Prw = 5.0 (R4B), and (h.l) 
PrM = 10.0 (R5B) for decaying MHD turbulence. 



Probability distribution functions of the local cross helicity He = u • b are shown 

The arguments of these PDFs are scaled by their 



via green full lines in Fig. 18 



standard deviations, namely, <Jhc'i data for the PDFs are obtained at tc for runs Rl- 
R5 in Figs. [T8](a.l)-[T8te.l), runs R3B-R5B in Figs. IT8yf.l)-[T8](h.l), and runs R1C-R4C 



in Figs. [I8[a.2)-|l8[d.2) for decaying MHD turbulence. For statistically steady MHD 
turbulence these PDFs are shown in Figs. 18 a.3)-[l8|^d.3) for runs R1D-R4D. All these 
PDFs have peaks close to He = 0; this reflects the tendency for u and b to be aligned 
or antialigned that we have discussed above. However, these PDFs are quite broad 
and distinctly non-Gaussian; this can be seen easily from the values of the mean fiHc^ 
standard deviation o"//^, skewness 'js^Hcy kurtosis 'J^^Hc given in Table |2| Thus 
fluctuations of He away from the mean are very significant. Table |2] also gives the value 
of the mean energy E and the ratio E/ hhci which does not appear to be universal; for 
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Figure 14. Semilog (base 10) plots of the PDFs cosines of angles, denoted generically 
by 9, between the current density j and the eigenvectors of fluid rate-of-strain tensor 
S, namely, (red full hue), (green dashed line), and (blue dotted line): (a.l) 
PrM = 0.1 (Rl), (b.l) PrM = 0.5 (R2), (c.l) PrM = 1.0 (R3), (d.l) Ptm = 5.0 (R4), 
(e.l) PrM = 10.0 (R5), (f.l) PrM = 1.0 (R3B), (g.l) PrM = 5.0 (R4B), and (h.l) 
PrM = 10.0 (R5B) for decaying MHD turbulence. 




cosS cosS cosO casff 



Figure 15. Semilog (base 10) plots of the PDFs of cosines of angles, denoted 
generically by 9, between the velocity u and the eigenvectors of the fluid rate-of-strain 
tensor S, namely, (red full line), (green dashed hne), and (blue dotted line): 
(a.l) PrM = 0.1 (Rl), (b.l) PrM = 0.5 (R2), (c.l) PrM = 1.0 (R3), (d.l) PrM = 5.0 
(R4), (e.l) PrM = 10.0 (R5), (f.l) PrM = 1.0 (R3B), (g.l) PrM = 5.0 (R4B), and (h.l) 
PrM = 10.0 (R5B) for decaying MHD turbulence. 

the runs R1-R5 and R3B-R5B it lies in the range 0.23-0.26, for R1C-R2C in the range 
—0.04-0.04, and for R1D-R4D in the range 0.05-0.2. For all our runs, with the exception 
of R2C, the mean fiHc the skewness 73,^^ cire positive. Even if the PDF of He had 
been a Gaussian, its mean value would have been within one standard deviation of 0; 
the actual PDF is much broader than a Gaussian. On symmetry grounds there is no 
reason for the system to display a nonzero value for fiHc unless there is some bias in 
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Figure 16. Semilog (base 10) plots of the PDFs of cosines of angles, denoted 
generically by 6*, between the magnetic field b and the eigenvectors of fluid rate- 
of-strain tensor S, namely, e\ (red full line), (green dashed line), and (blue 
dotted line): (a.l) PrM = 0.1 (Rl), (b.l) Prw = 0.5 (R2), (c.l) PrM = 1.0 (R3), (d.l) 
PrM = 5.0 (R4), (e.l) PrM = 10.0 (R5), (f.l) Prw = 1.0 (R3B), (g.l) PrM = 5.0 (R4B), 
and (h.l) PrM = 10.0 (R5B). 




Figure 17. Semilog (base 10) plots of PDFs of cosines of angles, denoted generically 
by 9, between (a) u and b, (b) u and w, (c) u and j, (d) us and j, (e) b and a;, and 
(f) b and j for runs Rl (red lines), R2 (green lines), R3 (blue lines), R4 (black lines), 
and R5 (cyan lines). 

the forcing or in the initial condition (the latter for the case of decaying turbulence). 
In any given run, if there is some residual He-, it is reflected in a slight asymmetry in 
alignment (or antialignment) of u and b, which we have studied above via the PDF of 
the cosine of the angle between u and b. When we consider the ratio hhc/^ it seems 
to be substantial in some runs but, given the arguments above, we expect it to vanish 
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Figure 18. Semilog (base 10) plots of the PDFs of the cross helicity Hq = u-b for (a.l) 
PrM = 0.1 (Rl), (b.l) PrM = 0.5 (R2), (c.l) PrM = 1.0 (R3), (d.l) Ptm = 5.0 (R4), 
(e.l) PrM = 10.0 (R5), (f.l) Prw = 1.0 (R3B), (g.l) Prw = 5.0 (R4B), (h.l) PrM = 10.0 
(R5B), (a.2) PrM = 0.01 (RIC), (b.2) PrM = 0.1 (R2C), (c.2) PrM = 1.0 (R3C), and 
(d.2) PrM — 10.0 (R4C) for decaying MHD turbulence; and for statistically steady 
MHD turbulence (a.3) PrM = 0.01 (RID), (b.3) PrM = 0.1 (R2D), (c.3) PrM = 1.0 
(R3D), and (d.3) PrM = 10.0 (R4D); the arguments of the PDFs are scaled by their 
standard deviations (THc- 



in runs with a very large number of collocation points; indeed, it is very small in runs 
R1C-R4C. 

Consider now the PDFs of the eigenvalues (blue dotted line), (green dashed 
line), and (red full line) of the rate-of-strain tensor S shown in Figs. [l9[a.l)- 



[T9[e.l) R1-R5 and Figs. [T9](f.l)-[T9](h.l) for runs R3B-R5B. Recall that these eigenvalues 
provide measures of the local stretching and compression of the fluid; also we label 
the eigenvalues such that A^ > A^ > A^. The incompressibility condition yields 
Yln=i = O5 whence it follows that A^ > and A^ < 0; the intermediate eigenvalue 
A^ can be either positive or negative. The illustrative plots in Figs. 19^a.l)-19|^h.l) 
from our decaying-MHD-turbulence runs show that the PDFs of A^ and A^ have long 
tails on the right- and left-hand sides, respectively. These tails shrink as we increase 
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Table 2. The mean i standard deviation i skewness ^^3,,Hc j ^^^d kurtosis fi4_Hc 
of the PDF of the cross hehcity He for our runs R1-R5 and R3B-R5B for decaying 
MHD turbulence at cascade completion; columns 6 and 7 give, respectively, the mean 
energy E and ratio of the means of the cross helicity and the energy, i.e., fiHc/^- 



r? n Tl 


1 1 TT 


(JHc 


/ / TT 


/ / J TT 


E 


llTJ IF, 


Rl 


0.118 


0.173 


1.103 


4.901 


0.461 


0.256 


R2 


0.118 


0.169 


1.096 


4.685 


0.467 


0.252 


R3 


0.120 


0.170 


1.096 


4.679 


0.490 


0.245 


R4 


0.112 


0.153 


1.003 


4.579 


0.477 


0.235 


R5 


0.105 


0.141 


0.934 


4.324 


0.460 


0.228 


R3B 


1.217 


1.804 
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4.912 
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R5B 


11.50 


17.05 


1.102 


5.000 


48.32 


0.238 


RIC 


0.014 


0.113 


0.615 


5.748 


0.358 


0.041 


R2C 


-0.224 


1.994 


-0.698 


8.441 


5.440 


-0.041 


R3C 


0.130 


2.005 


0.313 


5.637 


5.969 


0.022 


R4C 


0.859 


9.156 


0.364 


5.747 


29.05 


0.029 


RID 


0.169 


0.724 


0.737 


6.813 


3.090 


0.055 


R2D 


0.478 


0.886 


1.126 


5.954 


2.405 


0.199 


R3D 


0.454 


1.244 


1.904 


13.75 


3.039 


0.149 


R4D 


0.389 


1.110 


1.207 


9.225 


2.767 


0.140 
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Figure 19. Semilog (base 10) plots of PDFs of the eigenvalues (blue dotted line), 
(green dashed line), and A^ (red full line) of the rate-of-strain tensor § for (a.l) 
PrM = 0.1 (Rl), (b.l) PrM = 0.5 (R2), (c.l) Prw = 1.0 (R3), (d.l) PrM = 5.0 
(R4), (e.l) PrM = 10.0 (R5), (f.l) Prw = 1.0 (R3B), (g.l) PrM = 5.0 (R4B), and (h.l) 
PrM = 10.0 (R5B); the arguments of the PDFS are scaled by their standard deviations. 
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Pim [Figs. 19|^a.l)-(l9[e.l) for runs R1-R5, respectively], by increasing u while holding 
the initial energy fixed; thus, there is a substantial decrease in regions of large strain. 
However, if we compensate for the increase in u by increasing the energy in the initial 
condition such that k^a^rf^ and k^axVd both ~ 1, we see that these tails stretch out, 
i.e., regions of large strain reappear. 




Figure 20. Semilog (base 10) plots of PDFs of the local kinetic-energy dissipation rate 
eu (blue dashed line) and the magnetic-energy dissipation rate e& (red full line), with 
the arguments scaled by their standard deviations, for: (a.l) PrM — 0.1 (Rl), (b.l) 
PrM = 0.5 (R2), (c.l) PrM = 1.0 (R3), (d.l) PrM = 5.0 (R4), (c.l) PrM = 10.0 
(R5), (f.l) PrM = 1.0 (R3B), (g.l) PrM = 5.0 (R4B), (h.l) PrM = 10.0 (R5B), 
(a.2) PrM = 0.01 (RIC), (b.2) PrM = 0.1 (R2C), (c.2) PrM = 1.0 (R3C), and (d.2) 
PrM = 10.0 (R4C) for decaying MHD turbulence; and for statistically steady MHD 
turbulence (a.3) Ptm = 0.01 (RID), (b.3) PrM = 0.1 (R2D), (c.3) PrM = 1.0 (R3D), 
and (d.3) Ptm = 10.0 (R4D). 



We show PDFs of the kinetic-energy dissipation rate eu (blue dashed lines) and the 
magnetic-energy dissipation rate eb (red full lines) that are obtained at tc for runs R1-R5 
in Figs. |20Fa.l)-[20te.l), runs R3B-R5B in Figs. [20tf.l)-|20th.l), and runs R1C-R4C in 



Figs. |20[a.2)-[20Kd.2) for decaying MHD turbulence; and for statistically steady MHD 
turbulence they are shown in Figs. [20|(a.3)-|20|d.3) for runs R1D-R4D. All these PDFs 
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Table 3. The mean fi^^, standard deviation a^^, skewness 73, c„, and kurtosis 74. of 
the PDFs of the local fluid energy dissipation eu, and their analogues for e?,, for runs 
R1-R5, R3B-R5B, R1C-R4C, and R1D-R4D. 



Run 


















74,^6 


Rl 


0.0048 


0.0096 


6.382 


75.069 





0302 


0.0550 


7.611 


144.86 


R2 


0.0109 


0.0187 


6.053 


70.566 





0255 


0.0527 


8.182 


121.47 


R3 


0.0141 


0.0226 


5.450 


52.884 





0233 


0.0566 


10.46 


204.37 


R4 


0.0231 


0.0284 


4.042 


28.306 





0160 


0.0397 


7.955 


97.662 


R5 


0.0273 


0.0302 


3.684 


24.559 





0130 


0.0315 


6.682 


64.206 


R3B 


0.4165 


0.7345 


5.941 


70.070 





6440 


1.5881 


9.029 


147.71 


R4B 


6.7843 


10.898 


5.343 


55.676 


4 


4541 


13.377 


9.672 


155.71 


R5B 


21.164 


32.438 


5.353 


59.163 


9 


8332 


31.177 


9.621 


151.64 


RIC 


0.0031 


0.0076 


18.45 


1620.0 





0566 


0.0632 


3.340 


22.270 


R2C 


0.2354 


0.5177 


7.599 


112.92 


1 


5655 


3.5169 


13.99 


981.17 


R3C 


0.8349 


1.6375 


6.841 


105.85 


1 


3186 


3.5524 


10.41 


205.54 


R4C 


14.208 


22.900 


5.535 


66.496 


7 


1624 


24.974 


13.33 


406.21 


RID 


0.0448 


0.0630 


5.799 


89.678 





8087 


1.1004 


4.587 


40.328 


R2D 


0.0601 


0.0933 


5.180 


51.311 





4995 


0.9808 


8.488 


142.97 


R3D 


0.2886 


0.4233 


5.230 


55.366 





6389 


1.3120 


6.755 


80.154 


R4D 


0.4498 


0.5536 


5.055 


58.832 





5037 


1.1391 


8.077 


129.65 



have long tails; the tail of the PDF for ef, extends further than the tail of that for e„ for 
all except the smallest values of PrM [Figs. [20|^a.l), [20|(a.2), [20|^a.3) for runs Rl, RIC, 
and RID, respectively]. This indicates that large values of are more likely to appear 
than large values of and, given the long tails of these PDFs, suggests that, except at 
the smallest values of PrM we have used, we might obtain more marked intermittency for 
the magnetic field than for the velocity field. Furthermore, as we expect, the tail of the 
PDF of e„ is drawn in towards small values of e„ as we increase PrM [Figs. 20|^a.l)-20 e.l) 
for runs R1-R5, respectively] while holding rj and the initial energy fixed. However, if 
we compensate for the increase in u by increasing the initial energy so that kma^Vd 
kma.xVd both ~ 1, we see that the tails of the PDFs of ei, and eu get elongated as we 
increase PrM, e.g., in Figs. [20|^f . l)-[20|^h. 1) for runs R3B-R5B, respectively. The values of 
the mean fj,^^, standard deviation cTf^, skewness 73,e„, and kurtosis ■y^^^^ of the PDFs of 
the local fluid energy dissipation e„ are given for all our runs and their counterparts for eh 
are given in Table [3) From these values we see that the right tails of these distributions 
fall much more slowly than the tail of a Gaussian distribution. 

Similar trends emerge if we examine the PDFs of the moduli of the vorticity and 
the current density, u (blue dashed lines) and j (red full lines), respectively: These are 
presented at tc for runs R1-R5 in Figs. [2Tla.l)-l2T](e.l), runs R3B-R5B in Figs. [2T](f.l)- 
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Figure 21. Semilog (base 10) plots of PDFs of the moduli of the local vorticity 
(blue dashed lines) and the current density (red full lines), uj and j, respectively, with 
the arguments of the PDFs scaled by their standard deviations, for (a.l) PrM — 0.1 
(Rl), (b.l) PrM = 0.5 (R2), (c.l) Ptm = 1-0 (R3), (d.l) PrM = 5.0 (R4), (e.l) 
PrM = 10.0 (R5), (f.l) PrM = 1.0 (R3B), (g.l) PrM = 5.0 (R4B), (h.l) PrM = 10.0 
(R5B), (a.2) PrM = 0.01 (RIC), (b.2) PrM = 0.1 (R2C), (c.2) Prw = 1.0 (R3C), and 
(d.2) PrM — 10.0 (R4C) for decaying MHD turbulence; and for statistically steady 
MHD turbulence (a.3) PrM = 0.01 (RID), (b.3) PrM = 0.1 (R2D), (c.3) Prw = 1.0 
(R3D), and (d.3) PrM = 10.0 (R4D). 



|2Tfh.l), and runs R1C-R4C in Figs. |2T[a.2)-|2T](d.2) for decaying MHD turbulence; and 
for statistically steady MHD turbulence they are shown in Figs. 21 |^a.3)-|2l|^d.3) for runs 
R1D-R4D. The tail of the PDF for j extends further than the tail of that for for 
all except the smallest values of PrM [Figs. |2l|a.l), 21 a.2), and [2l|(a.3) for runs Rl, 



RIC, and RID, respectively], so large values of j are more likely than large values of 
u. Thus, given that these PDFs have long tails, it is reasonable to expect that, except 
at the smallest values of PrM we have used, intermittency for the magnetic field might 
be larger than that for the velocity field. Moreover, the tail of the PDF of u is drawn 
in towards small values of u as we increase PrM [Figs. [2l|(a.l)-|2l|^e.l) for runs R1-R5, 
respectively] while holding t] and the initial energy fixed; but if, while increasing z/, we 
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Table 4. The mecLii /i-i^, standard deviation cr^j, skewness ^nd kurtosis '74, cj 

of 

the PDFs of the modulus of the local vorticity w, and their analogues for j, for runs 
R1-R5, R3B-R5B, R1C-R4C, and R1D-R4D. 



Run 






73,0; 


74,0; 






73,i 


74,i 


Rl 


3.817 


3.138 


2 


174 


10.25 


3.109 


2.340 


2 


206 


11.42 


R2 


2.680 


1.938 


1 


961 


9.235 


2.794 


2.222 


2 


525 


13.84 


R3 


2.189 


1.512 


1 


871 


8.598 


2.606 


2.215 


2 


873 


17.16 


R4 


1.312 


0.766 


1 


534 


6.859 


2.121 


1.870 


2 


903 


15.74 


R5 


1.022 


0.567 


1 


363 


6.098 


1.906 


1.695 


2 


789 


13.95 


R3B 


11.50 


8.706 


1 


949 


8.908 


13.27 


12.06 


2 


791 


15.61 


R4B 


21.30 


14.98 


1 


828 


8.246 


32.58 


34.11 


3 


243 


19.05 


R5B 


26.93 


18.23 


1 


770 


8.049 


47.08 


51.90 


3 


360 


19.79 


RIC 


25.31 


23.10 


2 


265 


11.17 


21.00 


18.47 


2 


456 


14.56 


R2C 


15.70 


13.07 


2 


062 


9.831 


18.35 


17.95 


2 


899 


17.01 


R3C 


21.68 


15.49 


1 


746 


7.910 


38.85 


45.50 


3 


563 


23.00 


R4C 


2.945 


2.661 


3 


027 


26.37 


1.446 


0.861 


1 


451 


6.816 


RID 


12.14 


8.700 


2 


297 


14.16 


5.336 


3.456 


1 


702 


7.683 


R2D 


14.26 


9.802 


1 


745 


7.821 


12.58 


9.544 


2 


336 


12.64 


R3D 


10.07 


6.544 


1 


644 


7.463 


13.88 


11.24 


2 


320 


11.36 


R4D 


4.119 


2.349 


1 


388 


6.526 


11.92 


10.47 


2 


429 


12.46 



also increase the initial energy so that k^i^Vd ^^maxVd — 1, we see that the tails of 
the PDFs of j and cu get stretched out as we increase PrM, e.g., in Figs. 21 'f.l)-21 l^h.l) 
for runs R3B-R5B, respectively. The values of the mean /i^^, standard deviation a^^, 
skewness 73,^, and kurtosis 74^^ of the PDFs of the modulus of the local vorticity u for 
all our runs and their counterparts for j are given in Table |4j From these values we 
see that the right tails of these distributions fall much more slowly than the tail of a 
Gaussian distribution. 

We move now to PDFs of the local effective pressure (green full lines), which are 
shown at tc for runs R1-R5 in Figs. |22|a.l)-|22|(e.l) and runs R3B-R5B in Figs. |22|(f.l)- 
[22|^h.l) for decaying MHD turbulence; for statistically steady MHD turbulence they 
are shown in Figs. 22 a.2-d.2 for runs R1D-R4D. The values of the mean /ip, standard 
deviation cXp, skewness 73^^, and kurtosis 74^^ of the PDFs of the local effective pressure 
p are given for these runs in Table |5} These have mean fip = but are distinctly 
non-Gaussian as can be seen from the values of 73_p and 74^^. Pressure PDFs are 
negatively skewed in pure fluid turbulence as we have mentioned above; however, for 
MHD turbulence we find that the PDFs of the effective pressure p can be positively 
skewed, as in runs R1-R5, R3B-R5B, and run R4D, or negatively skewed, as in runs 
R1D-R3D; negative skewness seems to arise at low values of PrM- 
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Figure 22. Semilog (base 10) plots of PDFs of local effective pressure fluctuations 
(green full lines), with the arguments of the PDFs scaled by their standard deviations, 
for (a.l) PrM = 0.1 (Rl), (b.l) Prw = 0.5 (R2), (c.l) PrM = 1.0 (R3), (d.l) Ptm = 5.0 
(R4), (e.l) PrM = 10.0 (R5), (f.l) Ptm = 1.0 (R3B), (g.l) Ptm = 5.0 (R4B), (h.l) 
PrM = 10.0 (R5B), for decaying MHD turbulence; and for statistically steady MHD 
turbulence (a.2) PrM = 0.01 (RID), (b.2) Ptm = 0.1 (R2D), (c.2) Pym = 1.0 (R3D), 
and (d.2) PrM = 10.0 (R4D). 



The scale dependence of PDFs of velocity increments provides important clues 
about the nature of intermittency in fluid turbulence. To explore similar intermittency 
in MHD turbulence [72], we present data for the scale dependence of PDFs velocity 
and magnetic-field increments. As mentioned above, these increments are of the form 
5a||(x, /) = a(x + 1, t) — a(x, t)] ■ |, with a either u or b, / = |1| the length scale, and x 
an origin over which we can average to determine the dependence of the PDFs of ^ay on 
the scale /; for notational convenience, such velocity and magnetic-field increments are 
denoted by 6u{l) and 6b{l) in our plots. These PDFs are obtained at tc for runs R1-R5 
in Figs. |2ita.l)-[23te.l), runs R3B-R5B in Figs. [23tf.l)-|23th.l), and runs R1C-R4C in 



Figs. 23 a.2)-23|^d.2) for decaying MHD turbulence; and for statistically steady MHD 
turbulence they are shown in Figs. [23|^a.3)-|23|(d.3) for runs R1D-R4D. The PDFs of 
velocity increments are shown for separations / = 2Sx (red dashed thin line), / = lOSx 
(green dot-dashed thin line), and / = 1006x (blue full thin line), where 6x is our real- 
space lattice spacing; for PDFs of magnetic-field increments we also use the separations 
/ = 26x (black dashed line), / = 106x (cyan dot-dashed line), and / = 1006x (magenta 
full line); the arguments of these PDFs are scaled by their standard deviations. As in 
fluid turbulence, we see that these PDFs are nearly Gaussian if the length scale / is 
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Table 5. The mecin /z^, stciiidard. deviation fjp, skewness '73,p, and kurtosis "74. p of the 
PDFs of the local effective pressure p for runs R1-R5, R3B-R5B, and R1D-R4D. 



Run 






73,p 


74,p 


Rl 


-3.283E-16 


0.055 


0.224 


4.152 


R2 


-2.183E-16 


0.057 


0.256 


4.052 


R3 


6.606E-16 


0.061 


0.315 


3.722 


R4 


-2.787E-16 


0.060 


0.397 


3.493 


R5 


-6.596E-16 


0.059 


0.433 


3.527 


R3B 


2.975E-15 


0.609 


0.526 


5.283 


R4B 


1.323E-14 


3.184 


0.660 


5.645 


R5B 


-3.475E-14 


6.397 


0.719 


5.776 


RID 


9.014E-15 


0.738 


-0.533 


3.882 


R2D 


1.136E-15 


0.313 


-0.153 


4.697 


R3D 


1.244E-14 


0.589 


-1.066 


5.338 


R4D 


9.983E-16 


0.363 


0.221 


5.560 



large. As / decreases, the PDFs develop, long, non-Gaussian tails, a clear signature of 
intermittency. Furthermore, a comparison of the red and black dashed lines in these 
plots indicates that the PDFs of the magnetic-field increments are broader than their 
velocity counterparts in most of our runs; this suggests, as we had surmised from the 
PDFs of energy-dissipation rates given above, that the magnetic field displays stronger 
intermittency than the velocity field at all but the smallest values of PrM [Figs. [23|(a.l), 
[23|^a.2), and[23|^a.3) for runs Rl, RIC, and RID]; the general trend that we notice from 
these figures is that the magnetic-field intermittency is stronger than that of the velocity 
field at large magnetic Prandtl numbers but the difference between these intermittencies 
decreases as PrM is lowered. We will try to quantify this when we present structure 
functions in Subsection 13.51 

3. 5. Structure functions 

We continue our elucidation of intermittency in MHD turbulence by studying the scale 
dependence of order-p equal-time, velocity and magnetic-field longitudinal structure 
functions Sp{l) = (|5m||(x, and magnetic-field longitudinal structure functions 
Sp{l) = (|56||(x, respectively, where (5m||(x, /) = [u(x + l,t) — u(x, t)] ■ | and 
56||(x, /) = [b(x + — b(x, t)] ■ J. From these structure functions we also obtain 
the hyperfiatnesses F^{r) = S^{r)/[S^{r)f and F|(r) = S^{r)/[S^{r)f. For the inertial 
range r]^,ri^ -C Z -C L, we expect 3^(1) ~ l''p and Sp{l) ~ I'^p, where Cp and (p are 
inertial-range multiscaling exponents for velocity and magnetic fields, respectively; if 
these fields show multiscaling, we expect significant deviations from the K41 result 
qK4:1 _ (^bKAi = [Note that we do not expect any Iroshnikov-Kraichnan [73] scaling 
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iu.(;),«r(0 duUjMl) iu(l},6b{l] 

Figure 23. Semilog (base 10) plots of PDFs of velocity increments 6u{l), for 
separations I ~ 2Sx (red dashed thin line), 106x (green dot-dashed thin line), and 
lOOSx (blue, full thin line), and of magnetic-field increments Sb{l), for separations 
I = 2Sx (black dashed line), lOSx (cyan dot-dashed line), and lOOSx (magenta full 
line), with the arguments of the PDFs scaled by their standard deviations, for (a.l) 
PrM = 0.1 (Rl), (b.l) PrM = 0.5 (R2), (c.l) Ftm = 1-0 (R3), (d.l) Prw = 5.0 (R4), 
(e.l) PrM = 10.0 (R5), (f.l) Prw = 1.0 (R3B), (g.l) PrM = 5.0 (R4B), (h.l) PrM = 10.0 
(R5B), (a.2) PrM = 0.01 (RIC), (b.2) Pym = 0.1 (R2C), (c.2) Prw = 1.0 (R3C), and 
(d.2) PrM — 10.0 (R4C) for decaying MHD turbulence; and for statistically steady 
MHD turbulence (a.3) PrM = 0.01 (RID), (b.3) PrM = 0.1 (R2D), (c.3) PrM = 1.0 
(R3D), and (d.3) Prw = 10.0 (R4D). 

because we have no mean magnetic field in our simulations.] Given large inertial ranges, 
the multiscaling exponents can be extracted from slopes of log-log plots of structure 
functions versus /. However, in practical calculations inertial ranges are limited, so 
we use the extended-self-similarity (ESS) procedure [391 HQ] in which we determine the 
multiscaling exponent ratios QIQ cind Cp/Ch respectively, from slopes of log-log plots 
of (a) Sp versus 5*3 and (b) Sp versus 5*3; we refer to these as ESS plots. Our data for 
structure functions are averaged over 51 and 400 origins, respectively, for simulations 
with 512^ and 1024^ collocation points. 

We begin with data from our decaying-MHD-turbulence runs R1C-R4C, which use 
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Figure 24. Loglog (base 10) ESS plots of order-p structure functions of tiie 
velocity S^{1) [(a.l)-(d.l)] and magnetic-field S^{1) [(a.2)-(d.2)] versus S^{1) and 5^(0, 
respectively; plots of the local slopes of these curves are shown in the inset. The black 
horizontal lines, with vertical ticks at their ends, show the inertial range over which 
we have averaged the exponent ratios Cp/Cp and Cp/Cp! plots are shown for p = 1 
(red small-dotted line), p = 2 (green dot-dashed line), p = 3 (blue line), p = 4 (black 
thin-dashed line), p — 5 (cyan thick-dashed line), and p = 6 (magenta large-dotted 
line). Subplots (a.3)-(d.3) show the exponent ratios Cp/Ca versus p for the velocity 
(red dashed line with thin errorbars) and magnetic fields (blue dotted line with thick 
errorbars); the black solid line shows the K41 result Cjf^^^ = p/3. The semilog (base 10) 
plots (a.4)-(d.4) show the hyperflatnesses Fq{1) (red Hue) and Fq{1) (blue dashed line) 
versus I. Subplots in panels (a), (b), (c), and (d) are from our decaying-turbulence 
runs RIC, R2C, R3C, and R4C, respectively, with PrM = 0.01, 0.1, 1, and 10. 



1024^ collocation points and span the Pim range 0.01 — 10. Figures [24[a.l)-24|^d.l) 
show ESS plots for 5*^ (r) for runs R1C-R4C, respectively, for p = 1 (red small-dotted 
line), p = 2 (green dot-dashed hne), p = 3 (blue line), p = 4 (black thin-dashed line), 
p = 5 (cyan thick-dashed line), and p = 6 (magenta large-dotted line); their analogues 

h.2) 



for Sp{r) are given in Figs. 



24 



24 



d.2); the local slopes of these ESS curves are 
shown in the insets of these figures. Flat portions in these plots of local slopes help 
us to identify the inertial ranges. The regions that we have chosen for our fits are 
indicated by black horizontal lines with vertical ticks at their ends. In such a region, 
the mean value and the standard deviation of the local slope of the ESS plot for Sp{r) 
(or Sp{r)) yield, respectively, our estimates for the exponent ratio Q/Q {or Cp/Cs) 
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Table 6. Multiscaling exponent ratios Cp/Ca ^^'^ Cp/Ct from our decaying- MHD- 
turbulence runs R1C-R4C. 



p 


Q/Q; C/Cl(PrM = 0.01) 


Q/Q; C/Cl(PrM = 0.1) 


1 


0.41 ± 0.04 


0.35 ± 0.01 


0.39 ± 0.09 


0.42 ± 0.04 


2 


0.74 ± 0.04 


0.68 ± 0.01 


0.71 ± 0.08 


0.74 ± 0.03 


3 


1.00 ± 0.00 


1.00 ± 0.00 


1.00 ± 0.00 


1.00 ± 0.00 


4 


1.21 ± 0.09 


1.29 ± 0.03 


1.26 ± 0.14 


1.20 ± 0.03 


5 


1.38 ± 0.22 


1.56 ± 0.06 


1.51 ± 0.32 


1.37 ± 0.07 


6 


1.52 ± 0.41 


1.80 ± 0.10 


1.76 ± 0.53 


1.52 ± 0.13 


P 


Q/Q; C/Cl(PrM = 1) 


Q/Q; C/Cl(PrM = 10) 


1 


0.42 ± 0.03 


0.49 ± 0.04 


0.40 ± 0.06 


0.50 ± 0.05 


2 


0.74 ± 0.03 


0.80 ± 0.04 


0.72 ± 0.05 


0.81 ± 0.04 


3 


1.00 ± 0.00 


1.00 ± 0.00 


1.00 ± 0.00 


1.00 ± 0.00 


4 


1.25 ± 0.06 


1.15 ± 0.07 


1.27 ± 0.08 


1.15 ± 0.05 


5 


1.50 ± 0.16 


1.27 ± 0.18 


1.55 ± 0.18 


1.29 ± 0.12 


6 


1.74 ± 0.30 


1.38 ± 0.32 


1.83 ± 0.31 


1.45 ± 0.23 



and its errorbars. Figures [24[a.3)-|24[d.3) show plots of these exponent ratios versus 
p for the velocity field (blue dotted line with thick errorbars) and the magnetic field 
(red dashed line with thin error bars); the black solid line shows the K41 result for 
comparison. Though earlier studies [251 EH] have obtained such exponents from DNS 
studies, they have done so, to the best of our knowledge, only for PrM = 1; furthermore, 
they have not reported errorbars. Although our (conservative) errorbars are large, our 
plots of exponent ratios suggest the following: (a) deviations from the K41 result are 
significant, especially for p > 3, as in fluid turbulence; (b) at large values of PrM the 
magnetic field is more intermittent than the velocity field, in so far as the deviations of 
Cp/Cl from the K41 result p/3 are larger than those of Q/Q', (c) as we reduce PrM this 
difference in intermittency reduces until, at PrM = 0.01, the velocity field shows signs 
of becoming more intermittent than the magnetic field. This trend in intermittency is 
corroborated by plots versus I of the the hyperfiatnesses Fq{1) = ^^^^^3 (red line) and 

Fl^{l) = (blue dashed line) in Figs. 24 a.4)-24 ;d.4) for runs R1C-R4C, respectively: 
As I decreases, Fq{1) rises more rapidly than Fq{1 except at PrM = 0.01. 

Similar results follow from our studies of statistically steady MHD turbulence in 
runs R1D-R4D, which use 512^ collocation points and span the PrM range 0.01 — 10. 
Figures 25|^a.l)-25 d.l) show ESS plots for Sp{r) for runs R1D-R4D, respectively, for 
p = 1 (red small-dotted line), p = 2 (green dot-dashed line), p = 3 (blue line), p = 4 
(black thin-dashed line), p = 5 (cyan thick-dashed line), and p = 6 (magenta large- 
dotted line); their analogues for Sp{r) are given in Figs. 25 ^a.2)-25 d.2); the local slopes 
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Figure 25. Loglog (base 10) ESS plots of order-p structure functions of tiie 
velocity S^{1) [(a.l)-(d.l)] and magnetic-field S^{1) [(a.2)-(d.2)] versus S^{1) and 5^(0, 
respectively; plots of the local slopes of these curves are shown in the inset. The black 
horizontal lines, with vertical ticks at their ends, show the inertial range over which 
we have averaged the exponent ratios Cp/Cp and Cp /Cp ! plots are shown for p = 1 (red 
small-dotted line), p = 2 (green dot-dashed line), p = 3 (blue line), p — 4 (black thin- 
dashed line), p — h (cyan thick-dashed line), and p — Q (magenta large-dotted line). 
Subplots (a.3)-(d.3) show the exponent ratios Cp/Cs versus p for the velocity (red dashed 
line with thin errorbars) and magnetic fields (blue dotted line with thick errorbars) ; the 
black solid line shows the K41 result C,^'^^ ~ p/3. The semilog (base 10) plots (a.4)- 
(d.4) show the hyperflatnesses Fq{1) (red line) and Fq{1) (blue dashed line) versus 
I. Subplots in panels (a), (b), (c), and (d), are from our statistically steady MHD- 
turbulence runs RID, R2D, R3D, and R4D, respectively, with PrM = 0.01, 0.1, 1, and 
10. 



of these ESS curves are shown in the insets of these figures. We obtain estimates for the 
exponent ratio Cp/Ca ^"^^ Cp/Cs ^ind their errorbars as in Fig. 24, Figures 25 a.3)-25|^d.3) 
show plots of these exponent ratios versus p for the velocity field (blue dotted line with 
thick errorbars) and the magnetic field (red dashed line with thin errorbars); the black 
solid line shows the K41 result for comparison. Plots versus / of the hyperflatnesses Fq{1) 
(red line) and Fg(/) (blue dashed line) are given in Figs. 25 a.4)-25 d.4) for runs RID- 
R4D, respectively. All the trends here as exactly as in the decaying-MHD-turbulence 



plots in Fig. 24 



Tables |6] and [7] summarise, respectively, our results for multiscaling exponent ratios 
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Table 7. Multiscaling exponent ratios Cp/Cs ^-iid Cp/Ct from our statistically steady 
MHD-turbulence runs R1D-R4D. 



p 


Q/Q; C/Cl(PrM = 0.01) 


Q/Q; C/Cl(PrM = 0.1) 


1 


0.38 ± 0.04 


0.37 ± 0.01 


0.42 ± 0.03 


0.52 ± 0.11 


2 


0.72 ± 0.04 


0.70 ± 0.01 


0.74 ± 0.02 


0.83 ± 0.09 


3 


1.00 ± 0.00 


1.00 ± 0.00 


1.00 ± 0.00 


1.00 ± 0.00 


4 


1.23 ± 0.08 


1.26 ± 0.03 


1.20 ± 0.04 


1.12 ± 0.12 


5 


1.41 ± 0.19 


1.50 ± 0.07 


1.36 ± 0.11 


1.24 ± 0.28 


6 


1.55 ± 0.33 


1.72 ± 0.12 


1.49 ± 0.20 


1.39 ± 0.54 


P 


Q/Q; CpVCKPrM = 1) 


Q/Q; C/Cl(PrM = 10) 


1 


0.39 ± 0.04 


0.47 ± 0.02 


0.37 ± 0.02 


0.51 ± 0.06 


2 


0.73 ± 0.04 


0.79 ± 0.01 


0.71 ± 0.02 


0.82 ± 0.06 


3 


1.00 ± 0.00 


1.00 ± 0.00 


1.00 ± 0.00 


1.00 ± 0.00 


4 


1.20 ± 0.10 


1.13 ± 0.02 


1.25 ± 0.05 


1.11 ± 0.08 


5 


1.36 ± 0.30 


1.24 ± 0.06 


1.47 ± 0.13 


1.18 ± 0.16 


6 


1.46 ± 0.52 


1.33 ± 0.12 


1.65 ± 0.26 


1.25 ± 0.24 



for our decaying-MHD-turbulence runs R1C-R4C and our statistically steady MHD- 
turbulence runs R1D-R4D. The trends of these ratios with PrM have been discussed 
above. By comparing corresponding entries in the columns and rows of these tables, 
we see that exponent ratios from decaying and statistically steady MHD turbulence 
agree, given our (conservative) error bars. Thus, at least at this level of resolution 
and accuracy, we have strong universality of these exponent ratios, for a given value of 
Prjvi, in as much as the ratios from decaying-MHD turbulence agree with those from the 
statistically steady case. The dependence on PrM will be examined in Sec. |4} 



3.6. Isosurfaces 

As we have mentioned in our discussion of fluid turbulence, isosurface plots of quantities 
such as u, the modulus of the vorticity, give us a visual appreciation of small-scale 
structures in a turbulent flow; in fluid turbulence, iso-w surfaces are slender tubes if u 
is chosen to be well above its mean value [65l 171] . For the case of MHD turbulence it 
is natural to consider isosurface plots [71] of u, the modulus j of the current density, 
energy dissipation rates, and the effective pressure. 

Isosurfaces of co are shown at tc for runs R1-R5 in Figs. |26ta.l)-|26Fe.l), runs R3B- 



R5B in Figs. |26[f.l)-|26[h.l), and runs R1C-R4C in Figs. |26[a.2)-[26[d.2) for decaying 
MHD turbulence; and for statistically steady MHD turbulence they are shown in 
Figs. |26ta.3)-[26td.3) for runs R1D-R4D; these isosurfaces go through points at which 



the value of oo is two standard deviations above its mean value (for any given plot). 
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Table 8. A comparison of multiscaling exponent ratios Cp/Cs £^nd Cp/Cl from 
our statistically steady and decaying MHD simulations and from decaying MHD 
simulations by Mininni and Pouquet (Ref. '3U'), for PrM = 1- 



p 


Q Ref. 


m 




Q/Q Ref. m 


Q/Q (R3D) 


C^/Q (R3C) 


1 


0.30 


0.40 


0.39 ± 0.04 


0.42 ± 0.03 


2 


0.55 


0.74 


0.73 ± 0.04 


0.74 ± 0.03 


3 


0.74 


1.00 


1.00 ± 0.00 


1.00 ± 0.00 


4 


0.91 


1.22 


1.20 ± 0.10 


1.25 ± 0.06 


5 


1.04 


1.39 


1.36 ± 0.30 


1.50 ± 0.16 


6 


1.17 


1.56 


1.46 ± 0.52 


1.74 ± 0.30 


P 


C Ref. 




C/ci Ref. m 


C/Cl (R3D) 


C'p/Cl (R3C) 


1 


0.36 


0.43 


0.47 ± 0.02 


0.49 ± 0.04 


2 


0.63 


0.76 


0.79 ± 0.01 


0.80 ± 0.04 


3 


0.83 


1.00 


1.00 ± 0.00 


1.00 ± 0.00 


4 


0.97 


1.16 


1.13 ± 0.02 


1.15 ± 0.07 


5 


1.07 


1.28 


1.24 ± 0.06 


1.27 ± 0.18 


6 


1.14 


1.36 


1.33 ± 0.12 


1.38 ± 0.32 



For Prjvi = 1 it has been noted in several DNS studies that such isosurfaces are 
sheets [SI [251 EH ES| that there is a general tendency for such sheet formation 
in MHD turbulence; our results show that this tendency persists even when PrM 7^ 1. 
The number of high-intensity isosurfaces of to shrink as we increase PrM [Figs. |26ta.l)- 



26 e.l) for runs R1-R5, respectively], by increasing u while holding the initial energy 



fixed. However, if we compensate for the increase in u by increasing the energy in the 
initial condition such that k^axVd k^a.yjfd both ^ 1, we see that high-cu sheets 
reappear [Figs. [26|(f.l)-|26|^h.l) for runs R3B-R5B, respectively]. These trends are also 
visible in our high-resolution, decaying-MHD-turbulence runs R1C-R4C [Figs. [26|^a.2)- 
26{d.2)] and the statistically steady ones, namely, R1D-R4D [Figs.|2ga.3)-|26{d.3)]. One 
interesting point that has not been noticed before is that some tube-type structures 
appear along with the sheets at small values of PrM as can be seen by enlarging 
Fig. |26fa.3) for run RID. 

Similar features and trends appear in isosurfaces of j that are shown at t^. for 



runs R1-R5 in Figs. [27](a.l)-|27[e.l), runs R3B-R5B in Figs. [27[f.l)-|27[h.l), and runs 
R1C-R4C in Figs. |27[a.2)-|27(d.2) for decaying MHD turbulence; and for statistically 
steady MHD turbulence they are shown in Figs.|27|^a.3)-[27|^d.3) for runs R1D-R4D; these 
isosurfaces go through points at which the value of j is two standard deviations above its 
mean value (for any given plot). Again the dominant features in these isosurface plots 
are sheets; their number goes down as Prn increases with v while the initial energy 
is held constant; but if this energy is increased, the number of high- intensity sheets 
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a.1 




r., d.1 




Figure 26. Isosurfaces of the modulus ui of the vorticity: (a.l) PrM = 0.1 (Rl), (b.l) 
PrM = 0.5 (R2), (c.l) PrM = 1.0 (R3), (d.l) PrM = 5.0 (R4), (e.l) PrM = 10.0 
(R5), (f.l) PrM = 1.0 (R3B), (g.l) PrM = 5.0 (R4B), (h.l) PrM = 10.0 (R5B), 
(a.2) PrM = 0.01 (RIC), (b.2) PrM = 0.1 (R2C), (c.2) PrM = 1.0 (R3C), and (d.2) 
PrM = 10.0 (R4C) for decaying MHD turbulence; and for statistically steady MHD 
turbulence (a.3) PrM = 0.01 (RID), (b.3) PrM = 0.1 (R2D), (c.3) Ptm = 1-0 (R3D), 
and (d.3) PrM = 10.0 (R4D); these isosurfaces go through points at which the value 
of LJ is two standard deviations above its mean value (for any given plot). 



increase. 



Isosurfaces of e„ are shown at tc for runs R1-R5 in Figs. 28 'a.l)-28|^e.l), runs R3B- 
R5B in Figs. |28{f.l)-|28{h.l), and runs R1C-R4C in Figs. |28{a.2)-gd.2) for decaying 
MHD turbulence; and for statistically steady MHD turbulence they are shown in 
Figs. [28|^a.3)-|28|^d.3) for runs R1D-R4D; the isosurfaces go through points at which the 
value of eu is two standard deviations above its mean value (for any given plot). Similar 
isosurfaces of eb are shown at tc for runs R1-R5 in Figs. [29|^a.l)-29 ^e.l), runs R3B-R5B 
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Figure 27. Isosurfaces of the modulus j of the current density: (a.l) PrM — 0.1 
(Rl), (b.l) PrM = 0.5 (R2), (c.l) Ptm = 1.0 (R3), (d.l) Ptm = 5.0 (R4), (e.l) 
PrM = 10.0 (R5), (f.l) PrM = 1.0 (R3B), (g.l) PrM = 5.0 (R4B), (h.l) Ptm = 10.0 
(R5B), (a.2) PrM = 0.01 (RIC), (b.2) PrM = 0.1 (R2C), (c.2) PrM = 1.0 (R3C), and 
(d.2) PrM = 10.0 (R4C) for decaying MHD turbulence; and for statistically steady 
MHD turbulence (a.3) Pym = 0.01 (RID), (b.3) PrM = 0.1 (R2D), (c.3) Pim = 1-0 
(R3D), and (d.3) PrM = 10.0 (R4D); these isosurfaces go through points at which the 
value of j is two standard deviations above its mean value (for any given plot). 



in Figs. |29ff.l)-|29|(h.l), and runs R1C-R4C in Figs. |29|(a.2)-|29|d.2) for decaying MHD 
turbulence; and for statistically steady MHD turbulence they are shown in Figs. [29|^a.3)- 
[29|^d.3) for runs R1D-R4D; the isosurfaces go through points at which the value of 
is two standard deviations above its mean value (for any given plot). Here too the 
isosurfaces are sheets; they lie close to, but are not coincident with, isosurfaces of u and 
j; changes in PrM affect these isosurfaces much as they affect isosurfaces of u and j. 



Isosurfaces of p are shown at tc for runs R1-R5 in Figs. 30 a.l)-|30](e.l) and runs 
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Figure 28. Isosurfaces of the local fluid energy dissipation rate e„: (a.l) PrM = 0.1 
(Rl), (b.l) PrM = 0.5 (R2), (c.l) Ptm = 1.0 (R3), (d.l) PrM = 5.0 (R4), (e.l) 
PrM = 10.0 (R5), (f.l) PrM = 1.0 (R3B), (g.l) PrM = 5.0 (R4B), (h.l) PrM = 10.0 
(R5B), (a.2) PrM = 0.01 (RIC), (b.2) Ptm = 0.1 (R2C), (c.2) Prw = 1.0 (R3C), and 
(d.2) PrM = 10.0 (R4C) for decaying MHD turbulence; and for statistically steady 
MHD turbulence (a.3) PrM = 0.01 (RID), (b.3) PrM = 0.1 (R2D), (c.3) Pim = 1-0 
(R3D), and (d.3) PrM = 10.0 (R4D); these isosurfaces go through points at which the 
value of eu is two standard deviations above its mean value (for any given plot). 



R3B-R5B in Figs. 30 f.l)-30|^h.l) for decaying MHD turbulence; and for statistically 



steady MHD turbulence they are shown in Figs. [30[a.3)-|30](d.3) for runs R1D-R4D; the 
isosurfaces go through points at which the value of p is two standard deviations above 
its mean value (for any given plot). The general form of these isosurfaces is cloud- 
type, to borrow the term that has been used for isosurfaces of the pressure in fluid 
turbulence [60]. Here also changes in Ptm affect these isosurfaces much as they affect 
isosurfaces of u and j, in as much as high-intensity isosurfaces are suppressed as PrM 
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Figure 29. Isosurfaces of the local magnetic-energy dissipation rate e;,: (a.l) 
PrM = 0.1 (Rl), (b.l) PrM = 0.5 (R2), (c.l) Ftm = 1.0 (R3), (d.l) Ftm = 5.0 
(R4), (e.l) PrM = 10.0 (R5), (f.l) Pym = 1.0 (R3B), (g.l) Ptm = 5.0 (R4B), and 
(h.l) PrM = 10.0 (R5B). (a.2) PrM = 0.01 (RIC), (b.2) Ptm = 0.1 (R2C), (c.2) 
PrM = 1.0 (R3C), and (d.2) PrM = 10.0 (R4C) for decaying MHD turbulence; and for 
statistically steady MHD turbulence (a.3) PrM ^ 0.01 (RID), (b.3) PrM = 0.1 (R2D), 
(c.3) PrM = 1.0 (R3D), and (d.3) PrM = 10.0 (R4D); these isosurfaces go through 
points at which the value of eb is two standard deviations above its mean value (for 
any given plot). 

increases via an increase in u, unless this is compensated for by an increase in the initial 
energy (in the case of decaying MHD turbulence) or ReA- 
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Figure 30. Isosurfaces of the local effective pressure p: (a.l) PrM = 0.1 (Rl), (b.l) 
PrM = 0.5 (R2), (c.l) PrM = 1.0 (R3), (d.l) PrM = 5.0 (R4), (e.l) PrM = 10.0 
(R5), (f.l) PrM = 1.0 (R3B), (g.l) PrM = 5.0 (R4B), and (h.l) Pym = 10.0 (R5B) 
for decaying MHD turbulence; and for statistically steady MHD turbulencs (a. 2) 
PrM = 0.01 (RID), (b.2) PrM = 0.1 (R2D), (c.2) PrM = 1.0 (R3D), and (d.2) 
PrM — 10.0 (R4D); these isosurfaces go through points at which the value of p is 
two standard deviations above its mean value (for any given plot). 



3. 7. Joint probability distribution functions 



In this Subsection we present three sets of joint PDFs that have, to the best of our 
knowledge, not been used to characterise MHD turbulence. The first of these is a 
QR plot that is often used in studies of fluid turbulence as we have discussed in 



Subsections 2^ and 3A the next is a joint PDF of u and j; and the last is a joint 
PDF of and e^. 

We show QR plots, i.e., joint PDFs of Q and R, via filled contour plots; these are 



obtained at t^. for runs R1-R5 in Figs. 31 a.l)-|3l](e.l), runs R3B-R5B in Figs. [3l](f.l)- 



STj^h.l), and runs R1C-R4C in Figs. 31 a.2)-[3T](d.2) for decaying MHD turbulence; and 
for statistically steady MHD turbulence they are shown in Figs. [3l|(a.3)-|3l|^d.3) for runs 
R1D-R4D; the black curve in these plots is the zero-discriminant line D = ^R'^+Q^ = 0. 
These QR plots retain overall, aside from some distortions, the characteristic tear-drop 
structure familiar from fluid turbulence (see Subsection 3.1 and Fig. [5|. If we recall 
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Figure 31. QR plots, i.e., joint PDFs of Q and R shown as filled contour plots on a 
logarithmic scale for (a.l) Prw = 0.1 (Rl), (b.l) Pim = 0.5 (R2), (c.l) Pru = 1-0 (R3), 
(d.l) PrM = 5.0 (R4), (e.l) Ptm = 10.0 (R5), (f.l) PrM = 1.0 (R3B), (g.l) Ptm = 5.0 
(R4B), (h.l) PrM = 10.0 (R5B), (a.2) PrM = 0.01 (RIC), (b.2) PrM = 0.1 (R2C), 
(c.2) PrM = 1.0 (R3C), and (d.2) PrM = 10.0 (R4C) for decaying MHD turbulence; 
and for statistically steady MHD turbulence (a.3) PrM = 0.01 (RID), (b.3) PrM = 0.1 
(R2D), (c.3) PrM = 1-0 (R3D), and (d.3) PrM = 10.0 (R4D). The arguments Q and R 
of the QR plots are normalised by (w^) and (w^)'^/^, respectively. The black curve is 
the zero-discriminant line D = ^R"^ + = 0. 



our discussion of QR plots in Subsection |2.2| and we notice that, as we increase PrM 



[Figs. 31 a.l)-|3lKe.l) for runs R1-R5, respectively] while holding rj and the initial energy 
fixed, there is a general decrease in the probability of having large values of Q and R, 
i.e., regions of large strain or vorticity are suppressed; this corroborates what we have 
found from the PDFs and isosurfaces discussed above. However, if we compensate for 
the increase in z/ by increasing the initial energy, or Re\, so that fcmax^^ kuiaxVd 
both ~ 1, we see that Q and R can increase again. Note that when PrM is very small as 
in run RID [Fig. [3l|^a.3)], the tear-drop structure is very much like its fiuid-turbulence 
counterpart Fig. [5} which might well correlate with the appearance of some tube-type 
structures in the u isosurface in enlarged versions of Fig. [26|^a.3). 

We now consider joint PDFs of u and j that are obtained at tc for runs R1-R5 
in Figs. |32|a.l)-|32|e.l), runs R3B-R5B in Figs. |32|f.l)-|32|h.l), and runs R1C-R4C in 
Figs. 32 a.2)-32|^d.2) for decaying MHD turbulence; and for statistically steady MHD 




Figure 32. Joint PDFs of uj and j shown as filled contour plots on a logarithmic 
scale for (a.l) Ptm = 0.1 (Rl), (b.l) PrM = 0.5 (R2), (c.l) PrM = 1.0 (R3), (d.l) 
PrM = 5.0 (R4), (e.l) PrM = 10.0 (R5), (f.l) PrM = 1.0 (R3B), (g.l) PrM = 5.0 
(R4B), (h.l) PrM = 10.0 (R5B), (a.2) PrM = 0.01 (RIC), (b.2) PrM 0.1 (R2C), 
(c.2) PrM = 1.0 (R3C), and (d.2) PrM = 10.0 (R4C) for decaying MHD turbulence; 
and for statistically steady MHD turbulence (a.3) PrM = 0.01 (RID), (b.3) PrM = 0.1 
(R2D), (c.3) PrM = 1.0 (R3D), and (d.3) PrM = 10.0 (R4D). The arguments of the 
joint PDFs are normalised by their standard deviations. 



turbulence they are shown in Figs. [32Ka.3)-|32[d.3) for runs R1D-R4D. All these joint 
PDFs have long tails; as we move away from PrM = 1 they become more and more 
asymmetrical. Furthermore, as we expect, the tails of these PDFs are drawn in towards 
small values of u and j as we increase PrM [Figs. 32 a.l)-32 e.l) for runs R1-R5, 
respectively] while holding t] and the initial energy fixed. However, if we compensate 
for the increase in z/ by increasing the initial energy or Re^, so that kmnxVd k^axVd 
are both ~ 1, we see that the tails of the PDFs get elongated again. 

In the end we consider joint PDFs of e„ and eb that are obtained at tc for runs Rl- 
R5 in Figs. |33|(a.l)-|33{e.l), runs R3B-R5B in Figs. |33ff.l)-|33|(h.l), and runs R1C-R4C 
in Figs. 33|^a.2)-33 d.2) for decaying MHD turbulence; and for statistically steady MHD 
turbulence they are shown in Figs. [33|^a.3)-|33|(d.3) for runs R1D-R4D. The trends here 
are similar to the ones discussed in the previous paragraph. In particular, these joint 
PDFs have long tails; as we move away from PrM = 1 they become more and more 
asymmetrical; and the tails of these PDFs are drawn in towards small values of e„ and 




Figure 33. Joint PDFs of e„ and £{, shown as filled contour plots on a logarithmic 
scale for (a.l) Ptm = 0.1 (Rl), (b.l) PrM = 0.5 (R2), (c.l) Ptm = 1.0 (R3), (d.l) 
PrM = 5.0 (R4), (e.l) PrM = 10.0 (R5), (f.l) PrM = 1.0 (R3B), (g.l) PrM = 5.0 
(R4B), (h.l) PrM = 10.0 (R5B), (a.2) PrM = 0.01 (RIC), (b.2) PrM = 0.1 (R2C), 
(c.2) PrM = 1.0 (R3C), and (d.2) PrM = 10.0 (R4C) for decaying MHD turbulence; 
and for statistically steady MHD turbulence (a.3) PrM = 0.01 (RID), (b.3) PrM = 0.1 
(R2D), (c.3) PrM = 1.0 (R3D), and (d.3) PrM = 10.0 (R4D). The arguments of the 
joint PDFs are normalised by their standard deviations. 



eh as we increase Ftm [Figs. 33 a.l)-[33](e.l) for runs R1-R5, respectively] while holding 
7] and the initial energy fixed. But, if we make up for the increase in u by increasing 
the initial energy or Re^ so that k^axVd k^ajjfd both ~ 1, we see that the tails 
of the PDFs get elongated again. 



4. Discussions and Conclusion 



We have carried out an extensive study of the statistical properties of both decaying 
and statistically steady homogeneous, isotropic MHD turbulence. Our study, which 
has been designed specifically to study the systematics of the dependence of these 
properties on the magnetic Prandtl number PrM, uses a large number of statistical 
measures to characterise the statistical properties of both decaying and statistically 
steady MHD turbulence. Our study is restricted to incompressible MHD turbulence; 
we do not include a mean magnetic field as, e.g., in Refs. [33]; furthermore we do not 
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study Lagrangian properties considered, e.g., in Ref. [3l]. In our studies we obtain 
(a) various PDFs, such as those of the moduh of the vorticity and current density, 
the energy dissipation rates, of cosines of angles between various vectors, and scale- 
dependent velocity and magnetic-field increments, (b) spectra, e.g., those of the energy 
and the effective pressure, (c) velocity and magnetic-field structure functions that can 
be used to characterise intermittency, (d) isosurfaces of quantities such as the moduli of 
the vorticity and current, and (e) joint PDFs such as QR plots. The evolution of these 
properties with PrM has been described in detail in the previous Section. 

To the best of our knowledge, such a comprehensive study of the PrM— dependence 
of incompressible, homogeneous, isotropic MHD turbulence, both decaying and 
statistically steady, has not been attempted before. Studies that draw their inspiration 
from astrophysics often consider anisotropic flows [7S1 U7\ EHl EHl EHl Ell IH21 IM] , flows 
that are compressible [35l [M] , or flows that include a mean magnetic field [33l EH ESI EE] . 
Yet other studies concentrate on the alignment between various vectors such as u and 
b as, e.g., in Refs. [281 ESI [52]; some of these include a few, but not all, of the PDFs we 
have studied; and, typically, these studies are restricted to the case PrM = 1. Some of 
the spectra we study have been obtained in earher DNS studies but, typically, only for 
the case PrM = 1; a notable exception is Ref. [37], which examines the PrM— dependence 
of energy spectra but with a relatively low resolution. References [HI E21 EZ] have also 
considered some PrM— dependence but not for low PrM- Isosurfaces of the moduli of 
the vorticity and current density have been obtained earlier [2S1 EH IZl] for the case 
PrM = 1. The PrM dependence of these and other isosurfaces is presented here for the 
first time. The joint PDFs we have shown above have also not been investigated before. 

Here we wish to highlight, and examine in detail, the implications of our study 
for intermittency. Some earlier DNS studies, such as Refs. [30], had noted that, for 
the case PrM = 1, the magnetic field is more intermittent than the velocity field. This 
is why we have concentrated on velocity and magnetic-field structure functions. Our 
study confirms this finding, for the case PrM = 1- This can be seen clearly from 
the comparison of our exponent ratios, for PrM = 1, with those of the recent DNS 
of decaying-MHD-turbulence in Ref. [30] in Table [s] the errorbars that we quote for 
our exponent ratios have been calculated as described in the previous Section; we have 
obtained exponent ratios for Ref. [30] by digitising [HH] the data in their plot [Fig. 3 of 
Ref. [3^] of multiscaling exponents versus the order p (error bars are not given in their 
plot). Thus, at least given our errorbars, there is agreement between our exponent ratios, 
both for decaying and statistically steady MHD turbulence, and those of Ref. pO] for 
PrM = 1. We note in passing that the latter DNS is one of decaying MHD turbulence 
but with a very special initial condition, which allows an effective resolution greater 
than that we have obtained; however, the initial condition we use in our decaying- 
MHD-turbulence DNS is more generic than that of Ref. [30j- It is our expectation 
that nonuniversal effects, associated with different initial conditions [261 SHj) naight 
not affect multiscaling exponent ratios, except if we use nongeneric, power-law initial 
conditions [26] in which E{k) grows with k (at least until some large-/c cutoff). 
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Direct numerical simulations of decaying-MHD-turbulence, e.g., those of Refs. |25l 
|30] . often average data obtained from field configurations at different times that are 
close to the time at which the peak appears in plots of the energy-dissipation rate. This 
is a reasonable procedure, for PrM = 1, because the temporal evolution of the system is 
slow in the vicinity of this peak. We have not adopted this procedure here because, as 
we move away from PrM = 1, the cascade-completion peaks occur at different times in 
plots of e„ and as we have discussed in detail in earlier sections of this paper. 

Let us now turn to the PrM— dependence of the multiscaling exponent ratios shown 
in Tables g and |7| and in Figs|24ja.3)-|24(d.3) and |25{a.3)-|25|(d.3). Even though our 
error bars are large, given the conservative, local-slope error analysis we have described 
in the previous Section, a trend emerges: at large values of PrM the magnetic field is 
clearly more intermittent than the velocity field, in as much as the deviations of from 
the simple-scaling prediction are stronger than their counterparts for (p. However, the 
velocity field becomes more intermittent than the magnetic field as we lower PrM. Could 
this result, namely, the dependence of our multiscaling exponent ratios on PrM, be an 
artifact? We believe not. As we have discussed above, dissipation ranges in our spectra 
are adequately resolved; furthermore, we have determined exponent ratios from a a 
rather stringent local-slope analysis, which is rarely presented in earlier DNS studies of 
MHD turbulence. Ultimately, of course, this PrM dependence of multiscaling exponents 
in MHD turbulence must be tested in detail in very-high-resolution DNS studies of 
MHD turbulence; such studies should become possible with the next generation of 
supercomputers. 

It is useful to note at this stage that a recent experimental study of MHD turbulence 
in the solar wind [56j provides evidence for velocity fields that are more strongly 
intermittent than the magnetic field; this study does not give the value of PrM. However, 
their data for multiscaling exponents are qualitatively similar to those we obtain at 
low values of PrM. Furthermore, PDFs of He have also been obtained from solar- 
wind data [55]; these are similar to the PDFs we obtain for He- Of course, we must 
exercise caution in comparing results from DNS studies of homogeneous, isotropic, 
incompressible MHD turbulence with measurements on the solar wind where anisotropy 
and compressibility can be significant; and, for the solar wind, we might also have to 
consider kinetic effects that are not captured by the MHD equations. 

The last point we wish to address is the issue of strong universality of exponent 
ratios. In the fluid-turbulence context such strong universality [HI US] implies the 
equality of exponents (and, therefore, their ratios) determined from decaying-turbulence 
studies (say at the cascade-completion time) or from studies of statistically steady 
turbulence. Does such strong universality have an analogue in MHD turbulence? Our 
data, for any fixed value of PrM in Tables [6] and [7], are consistent with such strong 
universality of multiscaling exponent ratios in MHD turbulence; but, of course, our 
large errorbars imply that a definitive confirmation of such strong universality in MHD 
turbulence must await DNS studies that might become possible in the next generation 
of high-performance-computing facilities. 
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